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FOREWORD 


This  effort  ia  part  of  an  overall  program  whose  object  is  the  investigation, 
development  and  application  of  techniques  which  may  be  utilized  to  predict  elec¬ 
tronic  circuit  and  system  reliability,  both  initial  (time  zero)  and  time  dependent. 
The  investigation  of  prediction  techniques  and  tools  (i.  e.  ,  linearization,  Monte 
Carlo,  computer  poles  and  zeros  analysis,  regression  analysis)  conducted  by  the 
contractor  during  a  preliminary  phase  of  the  program  is  described  in  RADC-TR- 
61-299,  "Mathematical  Simulation  for  Reliability  Prediction". 

It  was  the  purpose  of  the  second  and  current  phase  to  develop  transfer  func¬ 
tions  for  selected  circuits  and  to  investigate  means  for  deriving  system  transfer 
functions  by  combining  the  transfer  functions  of  constituent  circuits.  In  follow  on, 
in-house  work  now  underway,  RADC  is  using  the  above  circuits  to  empirically 
evaluate  each  of  the  prediction  techniques  investigated  in  phase  one. 

In  order  to  assure  that  the  circuit  functions  chosen  would  be  appropriate  to 
typical  circuit  functions  encountered  in  Air  Force  ground  systems,  the  RADC 
preferred  functional  divisions  (Report  RADC-TR-59-243)  were  chosen  as  test 
vehicles.  Utilization  of  these  partic\ilar  circuit  functions  lead  to  the  additional 
advantage  that  the  transfer  functions  developed  and  the  predictions  of  performance 
characteristics  (initial  and  time  dependent)  will  lend  themselves  readily  to  form 
the  nucleus  of  a  library  of  preferred  circuit  transfer  functions  with  predicted  ini¬ 
tial  and  time  dependent  performance  characteristics. 
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ABSTRACT 


This  report  presents  the  results  of  the  second  phase  of  the  study  program 
for  the  development  of  techniques  for  predicting  the  reliability  of  electronic  sys¬ 
tems  from  statistical  information  about  the  performance  of  system  components. 

The  transfer  functions  which  have  been  developed  in  this  phase  are  mathematical 
models  of  the  actual  systems  to  be  evaluated.  They  are  used  to  determine  system 
performance  when  component  characteristics  vary  from  nominal  values  as  a  result 
of: 

(1)  Manufacturing  and  handling 

(2)  Degradation  due  to  age 

(3)  Internal  and  external  random  stresses 

This  report  concludes  that  the  simulation  techniques  used  are  the  most 
efficient  for  the  purpose  and  have  the  following  distinct  advantages. 

(1)  They  provide  a  means  of  determining  the  sensitivity  of  circuits  and 
systems  to  the  variations  noted  above.  When  circuits  and  systems 
become  complex,  then  the  simpler  conventional  methods  (e.  g.  ,  differ¬ 
entiation)  are  no  longer  feasible  for  determining  circuit  sensitivity. 

At  this  point,  mathematical  simulation  techniques  provide  the  most 
efficient  and  indeed  the  only  means  of  analyzing  and  evaluating 
system  performance. 

(2)  They  provide  a  means  of  determining  not  only  the  tolerance  limits  but 
also  the  shape  of  the  underlying  frequency  distribution  of  the  particular 
components  used  in  the  circuit  or  system  under  investigation. 

(3)  They  provide  an  unambiguous  definition  of  failure. 

(4)  They  provide  a  means  of  assessing  the  relative  merits  of  competing 
systems. 

(5)  They  allow  the  transfer  function  to  be  expressed  conveniently  in  matrix 
form  so  that  conv-entional  routines  available  to  modern  computers  can 
be  utilized.  Hence,  the  analytical  expressions  of  system  performance, 
necessary  even  for  the  most  rudimentary  statistical  design  techniques 
are  unnecessary.  (Such  analytical  expressions  are  laborious  to  form). 

(6)  They  can  be  used  to  optimize  a  system  for  a  given  cost. 

The  circuits  analyzed  in  this  report  will  not  only  serve  as  vehicles  to  sub¬ 
stantiate  the  findings  of  the  first  contractual  phase,  but  will  in  addition  serve  as  a 
nucleus  for  a  library  of  preferred  designs  with  defined  lifetime  characteristics. 
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SECTION  I 


INTRODUCTION 


1.  1  DISCUSSION  OF  BASIC  CONCEPTS 


1.  1.  1'  Transfer  Functions 

A  transfer  function  is  a  mathematical  function  in  an  equation  which  relates 
the  characteristics  of  system  performance  to  the  characteristics  of  component 
performance  and  to  system  inputs.  This  is  expressed  mathematically  by: 


z. 

1 


Ft(Xi,X2. 


yv'fr 


(1) 


where  the  represent  system  performance  characteristics,  the  x^^  represent 
component  characteristics,  the  y^  represent  inputs  into  the  system,  and  F^  is 
the  transfer  function. 

One  elegant  method  for  actually  solving  for  the  transfer  functions  describ¬ 
ing  a  particvilar  circuit  is  that  presented  in  [1  ,  2]  .  Here  one  writes  the  circuit 
equations  in  compact  matrix  form,  for  example: 


(2) 


The  n  by  n  square  matrix,  called  the  ''transformation  matrix",  consists  of 
elements  which  are  functions  of  values  of  the  various  system  components,  such  as 
resistance,  capacitance,  inductance,  etc. 
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The  individual  elements  of  the  transformation  matrix  are  functions  of  the 
nominal  values  of  the  component  parameters  and  are  derived  from  the  circuit 
equations  describing  the  particular  circuit;  examples  of  such  derivations  may  be 
found  in  Sections  III  and  IV.  The  b^  are  voltage  or  current  constants  also  obtain¬ 
ed  from  the  circuit  equations.  The  elements  of  the  third  matrix  are  the  indi¬ 
vidual  transfer  functions  describing  system  performance.  Following  Gabriel 
Kron,  we  may  say  that  the  transformation  matrix  dominates  the  whole  problem  of 
reliability  prediction.  Indeed,  once  this  matrix  is  known,  the  task  of  finding  the 
unknown  transfer  functions  can  be  made  entirely  automatic  -  a  mere  mechanical 
routine.  If  the  transfer  matrix  is  non-singular,  then  the  unknown  transfer  matrix 
[z]  can  be  found  by: 

[z]=[A]-^[B]  (3) 

Equation  (3)  thus  enables  one  to  find  the  z^  described  in  (1)  in  a  convenient  and 
compact  form. 

This  basic  procedure  of  finding  transfer  functions  can  be  used  to  perform 
any  one  of  three  different  types  of  circuit  analysis: 

1.  DC  Analysis:  The  performance  criteria  are  the  quiescent  currents 
and  voltages. 

2.  AC  Analysis:  The  performance  criteria  are  the  functional  relation¬ 
ships  between  amplitude,  phase  and  frequency. 

3.  Transient  Analysis:  The  performance  criteria  are  rise  and  decay 
time,  pulse  width,  overshoot,  etc.  Such  characteristics  are  called 
transient  responses. 

These  three  types  of  analysis  together  with  the  use  of  the  transfer  functions 
to  perform  them  can  be  applied  at  any  level  of  system  or  circuit  complexity. 
Specifically,  transfer  functions  can  be  applied  to  analyze  the  performance  of  an 
individual  circuit,  of  a  subsystem  composed  of  individual  circmt,  of  a  system 
composed  of  subsystems,  and  so  on  ad  infinitum. 

1.1.2  Monte  Carlo 

If  a  system  is  made  up  of  a  large  number  of  components  whose  character¬ 
istics  are  described  in  terms  of  statistical  measures,  e.  g. ,  distribution  functions, 
and  whose  characteristics  may  change  in  time  in  a  probabilistic  way,  then  the 
system  performance  criteria  may  also  be  determined  in  a  probabilistic  way  as 
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functions  of  the  distributions  of  component  characteristics.  Thus,  the  problem  of 
reliability  prediction  consists  in  determining  the  probability  distributions  of  meas¬ 
ures  of  component  performance.  The  functions  relating  system  performance  to 
component  characteristics  can  be  found  in  two  ways,  analytically  (see  Appendix  IV) 
or  synthetically. 

Given  the  cumulative  distribution  functions  of  various  characteristics  of  the 
components  of  a  system,  an  integral  equation  can  be  formed  which  relates  the 
corresponding  density  functions  to  cumulative  distribution  functions  of  character¬ 
istics  of  over-all  system  performance.  When  the  system  or  the  circiiit  is  rela¬ 
tively  simple,  i.  e. ,  when  the  underlying  density  functions  are  normal  and  the  com¬ 
ponents  are  either  all  in  series  or  all  in  parallel,  then  the  integral  eqiiation  can  be 
solved  by  analytic  means  (see  Appendix  IV).  However,  as  soon  as  the  circvdt  be¬ 
comes  in  any  way  more  complex,  or  if  the  underlying  density  functions  are  no 
longer  normal,  then  analytic  means  are  no  longer  practical  and  statistical  methods, 
of  solving  this  integral  equation  must  be  applied.  The  statistical  sampling  methods 
used  to  solve  such  integral  equations  are  called  Monte  Carlo  methods  or  mathe¬ 
matical  simulation  methods.  Examples  of  the  use  of  this  method  can  be  fo\md  in 
Section  lU  and  Section  IV  of  this  report,  as  well  as  in  [3] .  For  a  general  introduc¬ 
tion  to  Monte  Carlo  methods,  see  also  [4] . 

Before  concluding,  it  may  help  to  consider  an  example  of  the  application  of 
Monte  Carlo  methods  to  the  solution  of  an  area  integral.  Suppose  a  closed  curve 
to  be  traced  on  a  plane  surface  in  the  form  of  a  square,  as  shown  below.  On  this 
surface,  a  coordinate  system  is  superimposed.  A  source  of  pairs  of  random  num¬ 
bers  is  available,  any  pair  corresponding  to  a  point  on  the  square,  with  the  condi¬ 
tion  that  the  random  pairs  are  uniformly  distributed.  This  means  that  a  random 
point  may  fall  anywhere  on  the  surface  with  equal  likelihood.  If  these  random  pairs 
are  classified  into  two  groups  according  to  whether  they  fall  inside  or  outside  of 
the  closed  curve,  then  the  ratio  of  the  area  inside  the  curve  to  the  total  area  is 
approximated  by  the  ratio  of  the  points  falling  inside  the  curve  to  the  total  number 
of  random  points.  The  accuracy  of  this  approximation  can  be  increased  by  in¬ 
creasing  the  munber  of  random  pairs.  This  example  could  obviously  be  extended 
to  three  or  more  dimensions.  This  illustrates  the  connection  between  Monte  Carlo 
methods  and  integrals,  since  integrals  are  used  to  express  areas,  volumes,  etc. 
[11]. 
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Figure  1 


.  Statistical  Solution  of  Area  Integral 
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SECTION  II 


PROCEDURES  AND  COMPUTER  APPLICATIONS 

2.  1  INTRODUCTION 

In  this  section  a  description  is  to  be  found  of  the  actual  manner  in  which 
concepts  and  techniques  of  transfer  functions  and  Monte  Carlo  methods  are  applied 
in  a  reliability  analysis  of  a  circuit  or  a  system  of  circuits.  Next,  the  role  of  the 
computer  in  this  analysis  is  described  followed  by  the  flow  charts  used  in  the  com¬ 
puter  programs,  together  with  relevant  explanations  as  to  their  use.  Finally, 
there  is  a  formalized  description  of  these  procedures  suitable  for  emunerating 
specifications  in  contract  use. 

2.  2  GENERAL  PROCEDURE 

The  following  is  a  brief  description  of  the  procedure  for  performing  a  re¬ 
liability  prediction  as  outlined  in  the  block  diagram  in  Figure  2. 

Step  1:  Deterministic  Analysis 

A  set  of  performance  criteria  is  first  received  which  specifies  how  the 
system  i  s  to  perform.  Based  on  these  criteria,  a  circuit  is  then  selected,  pref¬ 
erably  from  a  catalog  of  standard  circuits  of  proved  performance  and  reliability. 
An  equivalent  circuit  is  then  formed  which  represents  the  schematic  and  which  is 
then  analyzed  so  as  to  yield  the  circuit  equations.  From  these  equations  the  trans¬ 
fer  function  is  then  formed.  In  many  cases  these  functions  can  be  formulated  in 
the  form  of  a  matrix.  This  matrix  is  set  up  with  the  specified  nominal  values  of 
the  circuit  and  solved  by  the  computer.  At  this  point,  if  it  is  impossible  or  in¬ 
convenient  to  find  an  equivalent  circuit;  or  to  set  up  the  circuit  equations;  or  to 
solve  the  matrices,  a  regression  analysis  may  be  alternatively  used  to  determine 
the  transfer  functions.  Substituting  the  nominal  values  of  the  component  charac¬ 
teristics  into  the  transfer  functions,  one  then  obtains  the  nominal  values  of  the 
system  performance  measures  represented  by  each  transfer  function.  The  nominal 
values  of  the  system  performance  measures  are  then  compared  with  the  nominal 
system  performance  criteria.  If  no  suitable  agreement  occurs,  then  a  new 
circuit  must  be  chosen  and  the  entire  procedure  repeated  until  suitable  agreement 
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Figure  2.  Procedure  for  Performing  Reliability  Prediction 


is  obtained.  At  this  point,  an  optional  experimental  check  of  the  calculated  nomin¬ 
al  values  may  be  undertaken.  The  computations  as  well  as  the  laboratory  experi¬ 
ments  must  be  repeated  until  computed  and  experimental  nominal  values  of  per¬ 
formance  measures  agree. 

Step  2:  Probabilistic  Analysis 

Once  it  has  been  ascertained  that  the  transfer  functions  are  correct,  the 
computations  correct,  and  the  laboratory  experiments  accurate,  so  that  all  nominal 
values  for  system  performance  agree,  then  the  next  stage  of  analysis  begins. 
Random  variations  in  the  component  characteristics  will  give  rise  to  random 
variations  in  the  system  performance  measures.  At  this  stage,  ramdom  numbers 
are  first  generated  according  to  the  rule  which  describes  the  random  variation  in 
the  component  characteristics,  i.  e. ,  the  generation  of  these  random  numbers 
simulates  the  random  variation  which  the  component  characteristics  are  subject  to. 
These  generated  random  numbers  are  then  substituted  into  the  trauisfer  functions, 
which  are  then  evaluated  and,  in  turn,  yield  sets  of  random  n\imbers  whose 
variations  can  be  described  by  a  cumulative  distribution  function  (see  Appendix  1). 
[5,  6]  This  cumulative  distribution  function  describes  the  resxilting  random  varia¬ 
tions  in  the  system  performance  measures.  This  process  can  be  repeated  at 
intervals  to  give  a  description  of  how  the  system  operates  in  time.  At  this  point, 
the  cumulative  distribution  functions  describing  the  system  performance  measures 
can  be  examined  to  see  if  they  are  in  accordance  with  the  specified  performance 
criteria.  If  they  are  not,  the  components  are  changed  and  the  Monte  Carlo  process 
described  above  is  repeated  until  the  performance  criteria  are  finally  satisfied. 

2.  3  COMPUTER  APPLICATIONS 

In  this  section  a  description  is  to  be  found  of  the  role  of  the  computer  in  the 
analysis  of  the  performance  of  circuits  and  systems,  and  in  the  efficient  treatment 
of  the  transfer  functions,  both  in  the  deterministic  phase  of  analysis  and  in  the 
probabilistic  phase  of  analysis  using  random  numbers.  Following  this  there  is  a 
generalized  flow  chart  which,  together  with  relevant  explanations,  describes  the 
computer  program  in  a  general  way  applicable  to  any  of  the  circviits  or  systems 
discussed.  With  slight  modifications,  it  can  be  rearranged  to  apply  to  any  of  the 
specific  circuits  or  systems  in  this  report. 


7 


2,  3.  1  Role  of  Computer 
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In  the  mathematical  simulation  procedures  the  computer  plays  an  indis¬ 
pensable  part.  In  the  probabilistic  stage  of  analysis  where  random  niunbers  are 
generated  and  substituted  in  the  transfer  functions  of  a  solution,  it  is  impractical 
to  do  the  work  manually.  The  time-saving  speed  of  the  computer  is  also  most 
desirable  in  the  solution  of  the  matrices  in  the  deterministic  stage.  The  speed 
and  convenience  of  using  the  computer  not  only  save  time  and  money,  but  also 
take  Monte  Carlo  methods  out  of  the  realm  of  theoretical  possibility  into  practical 
reality. 

A  further  advantage  in  the  use  of  computers  in  this  analysis  is  that  one 
obtains,  so  to  speak,  a  built-in  reference  library  which  completely  describes  the 
system  and  to  which  one  has  immediate  access.  A  documented  record  of  circuit 
performance  exists  on  cards  or  tape  and  information  retrieval  techniques  can  be 
used  to  ascertain  immediately  the  particular  characteristics  of  system  or  com¬ 
ponent  performance.  Given  a  particular  amplifier  placed  in  a  particular  configura¬ 
tion,  for  example,  one  can  immediately  find  the  gain,  bandwidth  and  associated 
distributions. 

2.  3.  2  Generalized  Flow  Chart  and  Explanation 

Figure  3  shows  a  generalized  flow  chart  which  (with  sviitable  modifications) 
can  be  applied  to  any  one  of  the  specific  circuits  or  systems  in  this  report.  The 
general  computer  routine  in  the  flow  chart  is  as  follows:  Initially,  the  machine  is 
prepared  for  operation  by  rewinding  the  magnetic  tape  and  setting  the  counter  to 
the  appropriate  number.  Next,  the  various  sets  of  random  numbers  are  generated 
in  accordance  with  the  underlying  frequency  distributions  of  the  component  param¬ 
eters.  This  involves  a  proper  adjustment  of  the  random  numbers,  i.  e. ,  fitting 
them  into  the  proper  scale,  e.  g. ,  0  to  1,  2.  5  to  4.  6,  900  to  1100,  and  ascertain¬ 
ing  whether  they  fall  within  the  proper  range  within  the  scale,  if  the  distribution 
is  to  be  truncated.  Next,  the  random  numbers  are  used  to  generate  and  solve  the 
matrix  involved  in  the  problem.  The  solutions  are  converted  for  the  plotter  and 
the  answers  written  on  the  magnetic  tape.  The  large  loop  on  the  left  indicates  that 
the  above  steps  are  repeated  one  at  a  time  for  each  single  random  number  and  that 
the  process  comes  to  a  halt  only  after  all  the  generated  random  numbers  have  been 
utilized  (counter  >  0).  [19] 
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Figure  3.  Generalized  Flow  Chart 


2.  4  FORMAL  PROCEDURE  FOR  RELIABILITY  ANALYSIS 


Step  1.  A  standard  circuit  or  a  system  composed  of  standard  circuits  is 
selected  to  satisfy  the  required  performance  criteria.  In  other  iwords,  the  circuit 
configuration  or  particular  arrangement  of  the  components  is  chosen,  as  opposed 
to  the  specific  or  nominal  values  and  their  respective  tolerances.  If  two  or  more 
competing  circuits  exist  and  are  capable  of  performing  similar  functions,  an 
analysis  may  be  made  of  them  in  order  to  determine  the  circuit  or  system  capable 
of  best  meeting  the  requirements. 

Step  2.  Next,  the  transfer  functions  describing  each  critical  performance 
criterion  such  as  bandwidth,  rise  time,  qviiescent  current  and  voltage,  noise 
figure,  etc. ,  are  generated.  The  transfer  functions  are  constructed  in  the 
following  manner: 

First,  equivalent  circuits  are  constructed  from  the 
schematics,  and  the  mathematical  equations  describing 
the  equivalent  circuits  are  formulated.  From  these 
circuit  equations,  the  matrices  contained  in  the  equation 
for  the  transfer  function  [Z]  =  [A]”^  [B]  can  then  be 
formed. 

In  cases  where  it  is  not  technically  or  economically  feasible  to  define~the 
transfer  function  analytically,  it  may  be  defined  experimentally  using  regression 
analysis  techniques. 

Step  3.  Next,  the  matrix  equation  obtained  in  Step  2  must  be  solved  by  the 
computer  using  the  specified  nominal  values  of  the  component  characteristics. 

Step  4.  Once  the  solution  is  obtained  mathematically,  an  experimental 
check  on  the  accuracy  of  the  transfer  function  may  be  made  in  the  laboratory  by 
constructing  the  actual  circuit,  taking  measurements  of  specified  performance 
criteria  and  comparing  these  experimental  values  with  the  values  obtained  analy¬ 
tically  from  the  transfer  functions. 

Note:  It  is  important  to  note  at  this  time  that  the  concept  of  failure  must  be 
clearly  defined  since  such  a  definition  is  essential  for  the  application  of  mathema¬ 
tical  simulation  procedures.  In  this  program,  the  following  definition  applies: 
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A  failure  is  a  cessation  of  ability  to  perform  a  specified 
fvmction  or  functions  within  previously  established  limits 
on  specific  performance  characteristics,  i.  e.  ,  referring 
to  the  cumulative  distribution  function  plotted  in  Figure 
4,  a  failure  is  any  unit  that  falls  in  the  shaded  area. 


UMIT  LIMIT 

Figure  4,  Cumulative  Distribution  Function  of  Critical  Performance  Limits - 

Used  to  Define  Failures 

Step  5.  Random  numbers  are  generated  according  to  the  rule  which  des¬ 
cribes  the  random  variation  in  component  characteristics,  i.  e.  ,  the  generation  of 
these  random  numbers  simulates  the  random  variation  which  the  component 
characteristics  are  subject  to.  These  variations  may  be  due  to  external  stresses, 
e.  g.  temperature,  vibration,  humidity,  etc.  ,  as  well  as  internal  stresses. 

Note:  Random  number  tables  for  frequently  occurring  distrubutions  which  repre¬ 
sent  a  large  class  of  those  elements  comprising  electronic  circuits  and  systems 
are  available  in  reference  [3]  ,  which  describes  the  generation  of  random  numbers 
having  the  following  distributions; 

1.  Uniform 

2.  Exponential 

3.  Weibull 

4.  Normal  (Gaussian) 

5.  Log  Normal 

6.  Poisson 

7.  Chi-Square  with  even  degrees  of  freedom 
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For  cases  where  it  is  not  known  if  any  of  the  above  sets  of  random  numbers 
describe  the  component,  empirical  cumulative  distribution  functions  can  be  gen¬ 
erated.  The  values  of  the  actual  components  to  be  used  in  the  system  are  meas¬ 
ured  and  then  ordered  to  form  cumulative  distribution  functions.  These  values  can 
then  be  used  directly  or  can  be  obtained  from  a  curve  that  has  been  fitted  to  these 
points. 


Step  6.  The  appropriate  random  numbers  are  then  substituted  into  the 
transfer  function,  which  is  then  evaluated  on  the  computer.  This  process  results 
in  cumulative  distribution  functions  describing  system  performance  measures[  17, 18] . 

Step  6a.  If  the  results  of  Step  6  are  not  satisfactory,  the  components  used 
in  Step  5  must  be  replaced  or  alternatively  the  basic  design  must  be  modified  so 
that  system  performance  is  less  dependent  on  component  tolerance  (e.  g,  ,  add 
feedback,  temperature  stabilization,  etc.  ). 

Step  7.  Repeat  steps  1  through  6  as  required. 
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SECTION  m 


CIRCUIT  APPLICATIONS 

3.  1  INTRODUCTION  TO  CIRCUIT  APPLICATIONS 

It  is  the  purpose  of  this  section  to  present  an  individual  circuit  analysis  of 
each  of  the  circuits  that  make  up  the  systems  represented  in  Figures  5  and  6.  The 
Timing  Network  is  composed  of  conventional  circuits  chosen  from  the  list  of 
preferred  circuits  included  in  the  RADC  Technical  Report  -  RADC-TR-59-243. 

The  Error  Sensing  and  Readout  System  is  composed  of  microminiature  circuits. 
Microminiature  circuits  were  chosen  since  they  are  of  a  "state  of  the  art"  nature 
and  consequently  relatively  little  is  known  of  their  reliability.  Therefore,  with 
the  cooperation  of  the  Sylvania  Microelectronics  Laboratory,  it  was  decided  to  in¬ 
clude  these  circuits  as  a  portion  of  the  Mathematical  Simulation  Program. 

It  is  possible  to  present  three  different  analyses  of  any  circuit;  namely  DC, 
AC  or  transient.  However,  since  the  purpose  of  this  program  is  to  prove  the 
feasibility  of  a  technique,  the  main  emphasis  has  been  placed  on  the  DC  or  quies¬ 
cent  current  and  voltage  portions.  The  techniques,  as  described  below  are  also 
applicable  to  the  AC  and  transient  analyses.  The  major  difference  in  the  approach 
is  due  to  the  fact  that  the  equations  derived  in  the  AC  or  transient  analyses  are 
time  dependent  or  differential  equations,  whereas  the  DC  analysis  yields  algebraic 
equations.  However,  it  is  possible  to  change  these  time  dependent  expressions 
into  algebraic  expressions  by  use  of  the  Laplace  Transformation  and  then,  once 
the  solution  is  obtained,  transform  it  back  into  the  time  domain.  [10] 

The  technique  which  has  been  referred  to  above  is  outlined  as  follows: 

Once  the  circuit  or  system  is  chosen,  a  complete  analysis  is  performed.  This  is 
realized,  in  the  DC  treatment,  by  describing  the  circviit  or  system  with  a  set  of 
nodal  and/or  mesh  equations.  These  equations  are  arranged  and  a  transformation 
matrix  is  formed.  Contained  in  this  matrix  are  all  of  the  critical  ciremt  param¬ 
eters  as  well  as  all  of  the  individual  components  and/or  various  combinations  of 
these  components  which  will,  if  varied,  cause  a  change  in  the  value  of  the  corres¬ 
ponding  parameter  or  parameters.  The  matrix  is  then  solved  for  the  nominal 
values  of  the  components  by  means  of  a  computer.  At  this  point  it  is  desirable  to 
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Figure  5,  Timing  Network 
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Figure  6.  Error  Sensing  and  Readout  System 


verify  the  accuracy  of  the  transfer  function  by  measuring  the  circuit  parameters 
on  a  breadboard  of  the  particular  circuit  or  system.  A  comparison  of  the  empiri¬ 
cal  and  analytical  values  will  accomplish  the  verification.  This  is  the  extent  to 
which  this  section  covers  this  technique;  however,  Section  IV  will  present  the  com¬ 
plete  solution,  the  remainder  of  which  is  briefly  described  below. 


Sets  of  random  numbers  are  generated  that  correspond  to  the  various  parts 
of  the  circuit  or  system.  The  type  of  random  numbers  depend  upon  the  type  of 
distribution  of  the  part  or  parts;  i.  e. ,  normal,  viniform,  etc.  These  random 
numbers  are  substituted  into  the  various  matrices  and  the  effects  on  the  per¬ 
formance  noted.  In  this  manner,  it  is  possible  to  isolate  the  most  critical  parts 
and,  therefore,  circuit  redesign  or  maintainability  steps  may  be  taken.  It  is  also 
possible  to  repeat  the  above  procedure  at  some  time  or  times  in  the  future  by 
generating  new  sets  of  random  numbers  based  upon  known  characteristics  of  the 
particular  parts,  thus  obtaining  a  prediction  of  circmt  reliability. 


3.  2  TRANSFER  FUNCTIONS  OF  CONVENTIONAL  CIRCUITS 

This  section  presents  the  transfer  functions  of  a  number  of  circuits  chosen 
from  the  list  of  preferred  circuits  included  in  the  RADC  Technical  Report  -  RADC- 
TR-59-243.  [9]  The  circuits  chosen  for  analyses  are  preferred  circuits  of  the 
highest  reliability  (failure  rate  wise)  which  have  undergone  extensive  life  test. 
Reference  Report  RADC-TR-59-243,  "Reliable  Preferred  Solid  State  Functional 
Divisions",  and  the  RADC  Reliability  Notebook  in  which  it  is  included.  In  this 
report  there  is  some  philosophy  about  the  circuits  which  you  may  find  helpful. 

3.  2.  1  Transfer  Function  of  the  Trigger  Circuit 

This  subsection  presents  an  analysis  of  an  emitter  coupled  binary  circuit, 
commonly  known  as  the  Schmitt  Trigger,  in  terms  of  its  transfer  functions.  [7] 

The  expressions  derived  are  those  for  the  circuit  appearing  on  page  47  of  the 
RADC  Technical  Report,  RADC-TR-59-243.  The  solutions  are  in  the  form  of 
matrices  of  algebraic  equations  for  the  quiescent  currents  and  voltages. 
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The  circuit  as  analyzed  is  shown  in  Figure  7  and  functions  as  a  general 
purpose  multivibrator  or  squaring  assembly.  The  transistors  Q2  and  per¬ 
form  as  emitter  followers  isolating  the  timing  circuit  from  external  loads;  and 
and  Q^,  coupled  through  the  common  emitter  resistor  perform  the 
squaring  fimction.  The  circuit  will  accept  a  sine  wave,  complex  wave,  or  rec¬ 
tangular  input  signal  and  present  two  DC  coupled  complementary  signals  at  the 
output  terminals. 

A  complete  solution  of  the  circuit  indicates  that  neither  transistor  when 
operating  was  in  a  saturated  condition.  A  verification  of  this  fact  may  be  re¬ 
ceived  by  an  examination  of  the  computer  solutions  of  the  base  currents  for  both 
cases  see  Tables  2and3.  From  the  solutions, it  may  be  seen  that  the  base  currents 
are  all  negative  indicating  an  impossible  situation.  Asa  check  on  the  accuracy  of 
the  transfer  functions,  a  breadboard  circuit  was  set  up  and  the  various  voltages 
were  measured.  As  may  be  seen  from  Table  1,  the  measured  and  calculated 
values  compere  very  closely,  thus,  confirming  the  accuracy  of  the  transfer 
functions. 


TABLE  1 

COMPARISON  OF  CALCUUTED  VS.  MEASURED  VALUES  OF  TRIGGER  CIRCUIT 

A 


NON-SATURATED 

,  Q,  OFF 

''e2 

''e3 

''e4 

CALCULATED 

22.1 

8.8 

7.4 

MEASURED 

22.5 

8.8 

8.0 

B 


Q,  NON-SATURATED 

,  OFF 

V 

''1 

''ei 

^^82 

''es 

CALCULATED 

10.0 

-0.7 

8.6 

8.7 

22.3 

MEASURED* 

11.0 

-0.7 

9.0 

8.6 

22.5 

*MMMirad  with  a  pMitIv*  DC  trigger  of  12  v. 
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Figure  7.  Trigger  Circuit  Schematic 
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Quiescent  Current  and  Voltage  for  Off 

When  the  condition  exists  that  Q,  is  off  and  Q.  is  in  either  the  saturated 
or  non-saturated  state,  ,  CR^  and  the  associated  circuitry  are  rennoved 
from  the  circuit  because  of  reverse  biasing.  Figure  8  is  the  equivalent  circuit 
used  for  the  circuit  analysis  and  Figures  9  and  10  are  the  equivalent  circuits  for  a 
transistor  in  the  non-saturated  and  saturated  state  respectively.  A  0.  6V  drop  was 
assumed  across  the  diodes  CR^  and  CR^  due  to  forward  biasing.  Utilizing  the 
nodes  indicated  in  Figure  8,  the  circuit  equations  were  derived  for  the  state  in 
which  is  saturated  and  are  presented  below. 
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The  matrix  used  for  the  computer  solution  of  this  state  is  presented  in 
equation  11.  The  computer  solution  of  this  matrix  is  included  in  Table  2.  For 
the  situation  where  is  non-saturated,  it  may  be  easily  seen  by  referring 


19 


A 


Figure  9.  Equivalent  Circuit  for 
Transiiton  in  Active  Region 


Figure  10.  Equivalent  Circuit  for 
Traniitton  In  Saturation 
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to  Figures  9  and  10  that  the  matrix  for  this  case  may  be  written  directly  from  the 
matrix  of  equation  1 1.  When  is  non- saturated,  the  base  current  is 

directly  proportional  to  l/P^  times  the  collector  current  Therefore,  the 

matrix  for  (non- saturated)  may  be  written  as  shown  in  equation  12.  The  com¬ 
puter  solutioix  for  this  matrix  is  shown  in  Table  3. 


TABLE  2 

COMPUTER  SOLUTION  OF  TRIGGER  CIRCUIT 
(Q^  OFF,  SATURATED) 


^12 

''e3 

''e4 

'C2 

'C3 

•C4 

'm  ' 

22.3 

8.3 

8.7 

0.004 

0.002 

0.002 

-0.0002 

TABLE  3 

COMPUTER  SOLUTION  OF  TRIGGER  CIRCUIT 
(Q^  OFF,  NON-SATURATED) 


^^82 

Ve4 

'C2 

•c3 

•C4 

22.1 

88 

24 

0.002 

0.001 

Quiescent  Current  and  Voltage  for  Off 


When  the  circuit  receives  a  sufficient  positive  DC  trigger  to  cause  to 
be  cut  off  and  to  conduct,  the  equivalent  circuit  represented  in  Figure  11 
applies.  In  this  case,  Q^,  CR^  and  CR^  are  ignored  due  to  reverse  biasing. 
The  first  situation  to  be  explored  is  when  is  saturated  and  the  equations 

for  this  case,  utilizing  the  circuit  indicated  in  Figure  11,  are  presented  in  equa¬ 
tions  13  through  20. 
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Figure  1 1 .  Equivalent  Circuit  for  off  and  Q,  on 


(18) 


V 


E3 


V 


E3 


24.  4 

(19) 

^C3  1 

^  P3  > 

(20) 

Its  corresponding  matrix  is  shown  in  equation  21.  The  computer  solutions 
for  this  matrix  for  values  of  V=1  to  10  volts  are  in  Table  4.  As  was  stated 
before,  it  is  possible  to  write  the  matrix  for  the  non- saturated  condition  directly. 
However,  in  this  case,  the  equations  are  as  shown  in  equations  22  through  28. 

The  matrix  for  this  set  of  equations  is  presented  in  equation  29  and  its  computer 
solution  in  Table  5. 
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(22) 


''ei  (iq)*  ‘01(13',)  = 
^E1  (•  157)  ^C1 

^(-i^)  "  '‘02(7*1 

''e3  (i-J  *  ‘C3  (^) 

“ei  (-i7  ■  ^(7)*  ‘ 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


The  analysis,  as  presented  above,  is  typical  of  the  analyses  of  the  remain¬ 
ing  circuits.  All  equations  necessary  for  the  complete  description  of  the  relation¬ 
ships  between  the  components  and  the  respective  circuit  performance  character¬ 
istics  are  defined  in  the  matrix  equations.  Thus,  all  required  tools  are  available 
for  a  Monte  Carlo  Mathematical  Simulation. 

3.  2.  2  Transfer  Functions  of  a  Variable  Frequency  Oscillator 

This  study  of  the  performance  of  a  variable  frequency  oscillator  will  be 
limited  to  an  analysis  of  the  oscillator's  frequency  as  a  fimction  of  its  components. 
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A  schematic  of  the  o&cxllator  is  shown  in  Figure  12  and  its  equivalent  circuit  is 
shown  in  Figure  13.  Table  6  prese-ts  the  performance  criteria  for  the  oscillator. 


TABLE  6L®J 

PERFORMANCE  CRITERIA  OF  VARIABLE  FREQUENCY  OSCILLATOR 


Requirements 

Maximum 

Design 

Center 

Minimum 

Frequency  range  t 

250  Icc 

____ 

50  kc 

Frequency  range  2 

1  me 

____ 

250  kc 

Frequency  range  3 

5  me 

1  me 

Frequency  range  4 

10  me 

5  me 

Power  Output 

4  mw 

3  mw 

2  mw 

Output  Impedance 

50 

— — 

Load  Impedance 

____ 

50  (2 

Output  Voltage,  50>ohm  Load  Open  Circuit 

Frequency  Stability,  Power  Supply 

l%v  DC 

— — 

DC  Supply  Voltage 

____ 

-12  V 

____ 

Operating  Temperature 

50OC 

25PC 

OOC 

It  has  been  shown  that  if  a  circuit  is  provided  with  a  sufficient  amount  of 
regenerative  feedback,  a  transistor  circuit  will  serve  as  a  generator  of  period¬ 
ically  varying  waves.  A  variety  of  feedback  circuits  which  differ  in  detail  are 
available  for  the  production  of  self- sustained  oscillations.  Characteristic  of  these 
circuits  is  a  feedback  network  through  which  a  portion  of  the  output  is  fed  back 
into  the  input  circuit  having  such  phase  and  amplitude  that  self-excitation  results. 
The  analysis  of  such  feedback  oscillators  will  be  simplified  with  the  application  of 
the  following  two  points:  [12,  15,  16,  25] 

1.  The  response  of  a  circuit  is  oscillatory  if  the  poles  of  the  transfer 
functions  are  complex  conjugates.  The  oscillation  will  be  sustained 
with  constant  amplitude  if  the  real  parts  of  these  poles  are  0. 

2.  The  voltage  amplification  or  transfer  function  of  a  feedback  amplifica¬ 
tion  with  voltage  feedback  is; 
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4>o  L-0.027h 

Cl-150/i/lF 
C2  *  1000/t/iF 
C3  *  1000/tflF 


NOTE:  THIS  CIRCUIT  IS  CONTAINED  IN  RADC  REPORT  TR59-243.  DEC.  15. 


Figure  12.  Variable  Frequency  Oscillator  Schematic 
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gm 

■«  ; 

;  R< 
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trc-aooo* 

C3>UQB^F  l*a027h 

WMi  THIS  CmCViT  IS  CONTAINIO  IN  RADC  UPORT  TRSf  20.  OCC.  », 


Figure  13.  Equivalent  Circuit  of  Variable  Frequency  Oscillator  at  50  kc 


0. 


A 
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The  poles  of  this  function  are  the  roots  of  the  denominator:  1  -  = 

The  procedure  for  analyzing  the  feedback  oscillator  is  as  follows; 

1.  Draw  the  complete  Class  A  equivalent  circuit  (this  is  shown  in  Figure 
13); 

a.  Replace  by-pass  capacitors  by  short  circuits 

b.  Replace  radio  frequency  chokes  with  open  circuits 

c.  Indicate  tube  or  transistor  capacitances 

d.  Draw  the  circuit  as  an  open  loop  amplifier 

2.  Compute  the  over -all  amplification  of  the  open  loop  circuit  as  a  function 
of  the  frequency  (p).  This  will  equal  the  feedback  factor;  A^. 

3.  Let  SA  -1  =  0. 

o 

4.  This  equation  will  appear  as  a  polynomial  in  (p)  with  unity  as  the 
coefficient  of  the  highest  powei  of  (p). 

5.  For  oscillations  to  occur  there  must  be  a  pair  of  conjugate  roots  with 
zero  real  parts.  These  two  conditions  yield  two  design  equations. 

From  Figure  13  a  nodal  analysis  allows  the  following  equations  to  be  written: 

V  -  V 

V,gl  t  (V,  -  +  -p— T  ''l  PS  =  “  '5»l 


(V, 


V,)g 


l'»b'e 


+  V 


284 


-  <^1- 


V,)g 


m 


=  0 


(31) 


(V  -  V  ) 

- r  ^  +  ^383  +  vs  =  0  <32) 


The  following  is  the  coefficient  matrix  of  the  above  equation: 
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pC- 


<8l'^8b'e'^  2  <-8b'e)  <“r 

p  LC,+  1  p  C 


-pC, 


p  ijCij+i  p  CjL+1 


Sb'e  -  <84  +  8b'e  + 


-pc, 

<-T-^ 


+  g_)  (-  g_)  (■ 


pC, 


p“LC,  +  1 


’m 


m 


p"'LC.  +  1 


+  g3  +  pCj) 


=  A 


The  gain  is  -  pA 
Vj  o 

V 

since  Vj  = 


(33) 

(34) 


where  is  the  coefficient  matrix  with  the  coefficients  of  replaced 
by  the  constant  terms  of  equations  (30,  31,  and  32)  respectively.  Then  the  gain, 

may  be  written  as: 

'^1 


V  V  ^ 
'^1  '^1 


(35) 


combining  equations  (33  and  35)  we  have 


“1 


-1  =  0 


multiplying  equation  (36)  by  V^A  and 

- 


VjA  =  0 


(36) 


(37) 


Since  A^  =  0  because  column  3  in  the  matrix  is  all  zeros  and  VI  cannot 
equal  zero,  then  A=  0.  (38) 

We  now  solve  the  coefficient  matrix  and  from  equation  (38)  we  set  the  matrix 
equal  to  zero.  Solving  the  matrix,  we  have 
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feb'eMs  *  8l*3«b'e  *  l*b' e*4S  *  S*b'e«4  *  '=2«b'.*3  * 

'=2«3*4  ■"  «ieb'e=l  ■"  ‘=3eieb’e  +  SiS^Cj  4  Cjgjg^  +  =3«b'e*3  *  '=3*3*4  * 

"=2*3  *  *l'=l  *  S*1  ■"  S*3  ^  '=3*4>*m'P  *  l‘=l'=2*b'=  +  '=2S*b'e  *  '=l'=2*4  * 


'=2'=3*4  '=l'=3*b'e  ■"  ‘=lS*4  "■  *b'e*4*3'=3‘-  *  *l*3*b'e‘=3'-  +  *1*3*4^''  * 


(39) 


(CjC,  4  CjCj  4  CjCj  4  gjgjCjWg  J  P  4  [g^.^g^CjC^L  4  gj^.^gjO^C,!.  4 
SC3g3g4L4gj*b,^C,C3L4gj*4C,C3L4(C,C3*3L4*jCjC3L)*J  4 


fSS'=3*b.,''"'=i'=2S*4'-*'=lSS'-  *ml'’  •“ 


The  circuit  values  for  this  matrix  are  shown  in  Table  1  and  substituting  these 
constants  in  the  above  equation  (39)  and  combining  terms  will  result  in  equation 
(40)  below. 


[0.  08457248618  x  lO"® 

+ 

2.  8  X  10"®  g  ] 

*m^ 

+[0.7996937653  x  lO"^^ 

+ 

3. 060708003  x  lO"^^  g  ] P 

•  20 

+[3. 5191363382  x  10  ^ 

+ 

4.769320769  x  10"^®  g  1 

m 

(40) 

+  [0.1233719939  x  lO”^^ 

+ 

0.1311538462  x  10"^^  g  ]  P® 

+  [0.1251125203  x  10’^°  +  0.405  x  10"^®  a  1  =  0 


The  first  technique  utilized  to  solve  equation  (40)  for  the  resonant  frequency 

(p)  was  graphical  in  nature.  A  computer  program  was  developed  to  solve  the 

above  equation  and  a  search  was  made  for  the  value  of  a  that  would  result  in  a 

®m 

pair  of  complex  conjugate  roots  on  the  imaginary  axis.  Figure  14  is  a  graph  of 
some  of  the  roots  of  equation  (40).  Figure  15  is  the  same  data  presented  on  an 
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expanded  scale.  This  graph  intersected  the  imaginary  axis  when  the  value  of 
was  0.  0086  and  p  »  j381,  300.  Since  p  *  jw  and  w*  2«£,  the  frequency  of  oscil¬ 
lation  may  be  found  by  equation  2jtfj  ■  381,  300j.  Divide  both  sides  by  j 


2irf  «  381,300 

£  =  ’  60,  717  CPS  (41) 

A  second  method  for  solving  equation  (40)  is  the  utilization  of  the  Routh-Hurwitz 

criterion.  The  technique  is  as  follows  for  a  4th  order  equation,  after  dividing 
•  8 

equation  (40)  by  10  the  4th  order  equation  appears  as  equation  (42) 

(.125x10"^^  +  .405  xl0^°g  )p^  +  (.123x10'^^  +  .  131  x  10"^^  g  )p^ 

m  ®m 

+  (3.52x10“^^  +  4.  77  X  10’^°  g^)p^  +  (.800  x  10'^  +  (42) 

3.  06  X  10'“*  g_)p  +  (.0846  +  2.  8  g  )  »  0 
“m  *m 

This  equation  is  of  the  same  form  as  equation  (43). 

1  -  A  B  «  ap^  +  bp^  +  cp^  +  dp  +  e  »  0  (43) 


The  formal  procedure  is  as  follows: 

1.  Form  the  array  from  the  coefficients 


A.  B.  .  , 

Ai^_  D.  -  m  .  B 

A  C 


(44) 


t 
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2.  The  number  of  sign  changes  in  the  first  column  of  the  array  equals  the 
number  of  roots  with  positive  real  parts.  If  there  are  no  sign  changes 
in  the  first  coliunn,  then  the  system  is  stable. 

3.  To  produce  a  sign  change  in  the  first  cobomn,  term  C  can  be  set  <  0. 

If  C  ■  0  is  solved,  it  yields  the  value  of  g  at  which  the  system 
becomes  luistable.  This  corresponds  to  the  condition  that  the  roots  of 
the  characteristic  equation  have  zero  real  parts. 

In  this  case: 

C  =  (.800  x  10"^+  3.  06x10"'*  g  )  (.  35  x  10"^®  +  .345  x  10"^® 

m 

g  -  .  62  X  10"^*  g  ^)  -  (8.  46  X  10"^  +  2.  8  g  )  (.  015  x  10"^^  + 

®m  *m  ®m 

.  0322  X  10"^°  g^  +  .  0172  X  10"^®  gj^^)  (45) 

This  can  be  solved  by  trying  positive  realizable  values  of  000  to  .  060) 

in  the  expression  to  find  the  s  at  which  C=0.  This  yields  a  value  of  g _  near 

“m  m 

.  007  or  .  008. 


(46) 


The  value  of  g_  at  which  C  -  0  can  be  substituted  into  the  characteristic 
m 

equation  and  this  can  be  solved  for  p. 

Since  we  are  interested  in  the  roots  having  zero  real  parts,  p  a  0  +  ju  ■  0 
+  j  2itf,  and  the  fourth  order  equation  can  be  reduced  to  two  equations: 

Re  =  0:  aw*  -  cu^  +  e  =  0  (47) 

Im  =  0:  -  bju^  +  dju  ■  0  (48) 

Both  equations  should  yield  the  same  value  for  f  when  the  proper  g^  is  used. 

The  frequency  of  oscillation  is  calculated  as  approximately  60  KC  as  shown  in 
Table  7  below. 


39 


TABLE  7 


COMPUTATION  OF  FREQUENCY  OF  OSCILUTION 


flm 

From  R«  =  0 

From  Im  =*  0 

0.007 

61.6  KC 

59  KC 

o.ooe 

61 .4  KC 

60  KC 

A  third  approach  relating  to  the  proceeding  results  is  an  explicit  determin¬ 
ation  of  the  resonant  frequency.  This  procedure  is  as  follows. 

Let  the  denominator  of  the  transfer  function  be  of  the  form  shown  in 
equation  (49)  where  the  coefficients  are  given  in  terms  of  g^  and  where  p  ■  ju 

A  +  bp  +  cp^  4  dp^  +  p^  =  0  (49) 


The  technique  which  follows  enables  one  to  solve  for  g  and  use  this 
result  to  find  the  frequency  of  oscillation.  Since  we  are  looking  for  the  complex 
conjugates,  it  is  possible  to  write  the  above  equation  with  these  factored  out. 
This  gives: 


or 


(P  +  iPg)  Ip  -  iPo>  (a  +  PP  +  P^)  =  0 
(P^  +  p^^)  (o  +  Pp  +  p^)  =  0 


Multiplying  out  and  equating  like  coefficients,  one  has 


(50) 


®Po^  +  (a  +  Po^)P^  PP^  P^  = 

a  +  bp  +  cp  +  dp  +  p 

(51) 

orPo^  =  a 

a  =  a/pQ^ 

(52) 

PPq^  =  b 

Pq^  =  ^/P 

(53) 

.  2 

o+Pq  =  c 

Of  +  b/p  a 

(54) 

P  *  d 

(55) 
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Combining  the  above,  we  have  the  following  two  equations: 


t  I  >  c  (56) 

Po"  ■  H  I”' 

From  these  two  results  the  solution  for  is  explicit  and  the  frequency 

of  oscillation  is  uniquely  determined. 

CONCLUSIONS: 

Of  the  three  methods  described  above,  it  appears  that  the  third  approach  is 
the  most  efficient.  The  results  determined  analytically,  satisfactorily  compared 
with  those  determined  empirically.  The  above  results  allow  a  Monte  Carlo  simu¬ 
lation  tD  be  performed  on  the  oscillator  circuit.  This  will  enable  the  distribution 
of  the  performance  parameters  to  be  determined  for  variations  in  circuit  compo¬ 
nents;  both  for  the  initial  performance  distribution  as  well  as  the  time -dependent 
distribution,  which  may  result  from  a  degradation  or  change  in  values  of  compo¬ 
nents. 

3.  2.  3  Transfer  Functions  for  the  High-Level  Amplifier 

•  The  next  circuit  to  be  analyzed  is  that  of  the  High-Level  Video  Amplifier 
shown  in  Figure  16.  This  amplifier  consists  of  a  two-stage,  common  emitter 
feedback  pair,  driving  a  common  emitter,  high  voltage  mesa  type  output  trans¬ 
istor.  Although  some  compression  may  occur  in  the  second  transistor,  the  entire 
chain  is  designed  as  a  linear  amplifier  up  to  the  output  capabilities  of  the  trans¬ 
istor  Q^,  Provision  is  made  for  DC  restoration  at  the  base  of  to  permit 
operation  with  long-duty  cycle  pulses.  However,  because  of  the  small  DC  offset 
inherent  in  this  form  of  clamp,  optimum  performance  is  achieved  only  if  the  duty 
cycle  is  restricted  to  20  percent  or  less. 

This  type  of  amplifier  is  designed  for  use  mainly  to  accept  positive  unipolar 
signals  from  low-level  sources  and  deliver  a  high-level  negative  control  voltage 
to  a  radar  or  some  type  of  cathode  ray  display  device. 
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Figure  16.  Video  High-Level  Amplifier  Schematic 


The  analysis  of  this  circuit  was  performed  in  two  parts;  the  first  being  that 
of  the  output  amplifier  or  Q^,  and  the  second  that  of  the  common  emitter  feed¬ 
back  pair. 


The  DC  equivalent  circuit  of  the  output  amplifier  is  shown  in  Figure 

17.  The  analysis  of  this  equivalent  circuit  yields  equations  (58) through  (&0)l 


^1  -  ^C3 


=  I 


C3 


^2-  Ve3--‘' 


^C3  ,  ^E3 

P3  ^10 


P3+  1 

'^E3  ■  ^C3  ^^12  ^13^ 


(58) 


(59) 


(60) 
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Figure  17.  DC  Equivalent  Circuit  of  Output  State  of  Video  High-Level  Amplifier 


The  equations  are  then  rearranged  and  put  into  matrix  form  for  computer 
solution.  When  typical  values  of  p  are  used,  the  matrix  equation  (61)  may  be 
solved  and  the  computer  solution  is  shown  in  Table  8. 


(61) 


COMPUTER  SOLUTION  OF  MATRIX  FOR  THE  OUTPUT  STAGE 


It  is  now  desirable  to  solve  for  the  DC  operating  point  for  the  common 
emitter  feedback  pair  (Q^  and  02)*  The  equivalent  circuit  for  this  problem  is 
shown  in  Figure  18. 


Figure  18.  Equivalent  Circuit  for  Common  Emitter  Pair  of  High-Level  Video  Amplifier 

The  circuit  analysis  of  this  equivalent  circuit  yields  equations  (62)  through 

(68). 
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Froii,  il.e  above  equations,  the  matrix  equation  (69)  is  lorined  and  the 
computer  solution  of  the  equation  is  presented  in  Table  9. 

TABLE  9 

COMPUTER  SOLUTION  OF  COMMON  EMITTER  PAIR  OF  HIGH-LEVEL  AMPLIFIER 


^Cl 

^C2 

Vei 

^E2 

^2 

'ci 

‘C2 

5.76 

16.78 

1.11 

5.16 

4.05 

0.003 

0.005 

Measurements  of  the  critical  circuit  parameter  were  made  on  a  breadboard 
of  the  amplifier  to  validate  the  transfer  functions  represented  in  equations  (61)  and 
(69).  A  comparison  of  the  analytical  and  empirical  results  is  presented  in  Table 
10. 


TABLE  10 

COMPARISON  OF  MEASURED  AND  COMPUTED  VALUES 
OF  VIDEO  HIGH-LEVEL  AMPLIFIER 


VOLTAGE 

COMPUTED 

MEASURED 

Vci 

5.756  V 

7.0  V 

^C2 

16.776 

14.5 

^El 

1.104 

0.95 

''e2 

5.156 

6.2 

^2 

4.055 

5.1 

97.78 

98.0 

0.313 

0.4 

As  may  be  seen  from  Table  10,  the  results  check  rather  closely.  The 
agreement  between  the  computed  and  measured  values  would  improve  if  the  meas¬ 
ured  values  of  and  (35  and  42  respectively)  were  used  in  place  of  the 
assumed  values  of  55  and  60.  However,  the  results  serve  to  validate  the  transfer 
functions  for  the  DC  bias  conditions. 


45 


Another  important  parameter  in  addition  to  the  DC  quiescent  point  is  that 
of  the  AC  small  signal  gain.  Therefore,  an  AC  equivalent  circuit  was  developed 
and  is  shown  in  Figure  19.  Since  the  frequency  AC  Equivalent  Circuit  of  Video 
High-Level  Amplifier  of  interest  is  50  KC,  the  circuit  may  be  simplified  to  that 
shown  in  Figure  20.  Utilizing  the  simplified  equivalent  circuit,  a  nodal  analysis 
yields  equations  (70)  through  (77).  These  equations  are  then  formed  into  matrix 
equation  (78).  This  matrix  is  then  solved  by  means  of  a  computer  and  the  results 
are  presented  in  Table  11. 

(V^  -  ^2)81  +  +  (V4  -  V^lgej  =  0  (70) 


P  l^bj  +  +  87)  -  P2^b2  =  ° 

-  V«5  *  \  (^^3)  =  » 

^484+  <^4  '  '^2)8®!  =  0 


(71) 

(72) 

(73) 


P3^b3  ^6811  "  ® 

ib^  =  (Vj  -V3)gj 


'k  "  Pl^b/  V3  = 


^b,  P3S, 


V8i2  ^  8^3/ 


Note:  g  =  ^ 


81 

.  001 

86,9,10  =  *00205 

811  ■ 

.  000454 

82 

= 

.0000213 

gej  =  .077 

Pi  = 

55 

83 

s 

. 000147 

gej  =  .077 

II 

CM 

CQ. 

60 

84 

sr* 

. 00555 

g^  *  . 00454 

P3  = 

13 

85 

s 

, 000213 

gl2  *  .0147 

(74) 

(75) 

(76) 

(77) 
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A-J.50  KC; 


TABLE  11 


COMPUTER  SOLUTION  OF  AC  CIRCUIT  ANALYSIS  OF 
VIDEO  HIGH-LEVEL  AMPLIFIER 


^^2 

^3 

m 

''5 

B 

'bi 

*b2 

*b3 

0.96 

|QQ|[ 

0.93 

14.6 

-379.8 

0.00004 

-0.0008 

0.008 

The  above  analysis  is  based  on  small  signal  equivalent  circuits  and  as¬ 
sumes  that  the  transistors  are  linear  devices.  In  practice  the  predicted  perform¬ 
ance  will  only  be  valid  for  low-level  signals.  Specifically,  the  output  stage 
will  be  extremely  non-lincar  to  large  input  voltages. 

Table  12  presents  the  computed  and  measured  values  of  small  signal 
voltage  gain,  input  impedance  and  output  impedance  in  the  50  KC  frequency  range. 

TABLE  12 

COMPARISON  OF  COMPUTED  AND  MEASURED  SMALL  SIGNAL  GAIN 
OF  VIDEO  HIGH-LEVEL  AMPLIFIER 


Transfer  Function 

or  Impedance 

Computed 

Measured 

Voltage  Gain 

V. 

V^=-380 
’  V 

-330  +  25% 

Input  Impedance 

18K 

Output  Impedance 

^out  *^11 

2.2K 

The  pulse  voltage  gain  of  the  amplifier  was  also  measured  using  a  10%  duty 
cycle.  This  transfer  function  is  plotted  in  Figure  21  and  illustrates  that  the  gain 
is  qviite  linear  for  less  than  100  MV  and  approximates  the  computed  small 

signal  value  in  this  region.  In  fact  from  Figure  21,  the  low-level  pulse  voltage 
gain  ■  ^cs^^n  *  38V/0.  IV  a  3B0.  This  indicates  that  the  small  signhl 
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voltage  gain  transfer  function  is  also  valid  for  the  low>level  pulse  voltage  gain  for 

duty  cycles  under  50%.  Transistor  begins  to  compress  the  signal  for  inputs 
above  100  MV  because  its  gain  drops  as  it  is  driven  to  cutoff. 

Table  13  presents  the  performance  criteria  of  the  amplifier. 

NOTE;  The  current  gain  and  the  base  emitter  voltage  of  the  three  transistors  used 
in  this  amplifier  were  taken  from  manufacturer's  documents;  [13]  and  [14]. 

3.  2.  4  Transfer  Function  of  the  Video  Low-Level  Limiting  Amplifier 

This  subsection  will  develop  the  transfer  function  for  the  DC  bias  conditions 
of  the  Video  Low-Level  Limiting  Amplifier  presented  on  page  3  of  the  RADC  Tech¬ 
nical  Report  RADC-TR- 59-243. 

The  circuit  accepts  short  duty  cycle  positive  pulses  from  a  radar  detector 
or  some  other  source  and  provides  a  nominal  voltage  gain  of  12  up  to  the  limiting 
level.  Input  signals  in  excess  of  5  volts  are  delivered  to  the  output  terminals  at  a 
5  volt  level.  The  maximum  signal  swing  capability  and  limiting  action  is  achieved 
by  biasing  the  first  stage  near  cutoff,  the  second  stage  near  saturation  and  the  out¬ 
put  stage  at  cutoff.  A  biasing  scheme  of  this  type  admits  unipolar  positive  pulses 
at  the  input  and  provides  for  minimum  standby  power  dissipation  in  the  amplifier. 
Thus,  power  dissipation  capabilities  of  the  output  stage  are  used  most  advantageous¬ 
ly  for  efficient  signal  power  transfer  to  succeeding  stages.  By  properly  control¬ 
ling  the  gain  of  the  first  two  stages,  the  output  stage  may  be  driven  to  saturation 
and  limiting  action  occurs  at  the  desired  input  signal  level.  The  performance 
criteria  of  the  amplifier  is  presented  in  Table  14. 

The  amplifier  schematic  is  shown  in  Figure  22.  Critical  performance 
criteria  that  are  considered  in  the  following  analysis  are  the  DC  quiescent  currents 
and  voltages  for  the  condition  of  no  input  signal. 

The  DC  equivalent  circuit  of  the  amplifier  is  presented  in  Figure  23.  Note 
that  all  AC  dependent  parameters  are  removed  from  the  circuit  and  the  diode  CRl 
is  removed  because  of  reversed  biasing.  The  analysis  of  this  circuit  yields 
equations  (79)  through  (87). 
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TABLE  14 


PERFORMANCE  CRITERIA  OF  VIDEO  LOW-LEVEL  LIMITING  AMPLIFIER 


REQUIREMENTS 

MAXIMUM 

DESIGN 

CENTER 

MINIMUM 

INPUT  LEVEL.  PEAK  POSITIVE 

5V 

— 

O.OOIV 

INPUT  IMPEDANCE 

— 

— 

1000  a 

OUTPUT  LEVEL  (75-OHM  LEAD) 

5V 

0.012V 

VOLTAGE  GAIN  (0.001V  TOO.  IV) 

15 

12 

8 

RISE  TIME 

0.05/iSEC 

0.035/1  SEC 

- - 

DROOP  (500-  SEC  PULSE) 

10% 

— 

__  —  _ 

OVERSHOOT 

5% 

- - 

- - 

LINEARITY  (0.001V  TO 0.05V) 

5% 

— 

_ - 

D-C  SUPPLY  VOLTAGE 

— 

+25V 

- - 

OPERATING  TEMPERATURE 

i _ 

85°C 

+12V 

-55°C 

'’eI  ,  .  ^C1  .  0.6  ,  ''^CZ  -  ^E1 
—  =  'ci  +  T7  "  “  - 


*'3 
25  -  V 


'1  *bl 


1C  Cl  p  r^2 


E2  ^E2  ■  ^E1  ■  ^  _  T  j.  ^C2  j.  0. 6 

6  - ^CZ  TT  "  ^ 


25-Vc2  ^  .  Vc2-Vej 

— ic —  =  ^C2  +  — 1C — 


(79) 

(80) 

(81) 

(82) 
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5S 


(83) 


Ve2  ^  0*^  =  ^Cl 


^C1  ^  _  ^E2  ~  ^E1  '  ^ 

^1  ^’bl 


^E3 

ITT 


I  +  i21 
C3  p- 


I 


C3 

P3 


12  -  Vgj  -  0.  6 


Ve3  +  0.  6 

— r: — 


12  -  V 


C3 


"C3 


10 


(84) 

(85) 

(86) 

(87) 


These  equations  are  used  to  form  the  matrix  equation  (88)> 

The  matrix  equation  is  then  oolved  using  a  computer  and  Table  15  presents 
a  comparison  of  the  measured  and  calculated  values  of  the  critical  circuit  param¬ 
eters.  Thus,  an  examination  of  the  resvilts  validates  the  developed  transfer 
function. 
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TABLE  15 

COMPARISON  OF  MEASURED  AND  CALCULATED  VALUES  OF  THE 
LOW-LEVEL  LIMITING  AMPLIFIER 


COMPUTED: 

''ci 

''a 

''e2 

‘ci 

*C2 

*C3 

11.76 

11.34 

11.68 

0.5656 

11. 16 

0.453 

0.000976 

0.000719 

MEASURED: 

''ci 

''C2 

^C3 

''ei 

'^E2 

H9 

11.0V 

12.0 

11.7 

0.55 

10.0 

HB 
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3.  2.  4.  1  Low-Level  Amplifier  Mathematical  Simulation 

As  was  described  in  detail  in  Section  Z,  the  performance  of  the  Low-Level 
Limiting  Amplifier  was  evaluated  by  simulating  the  assembling  of  1000  amplifiers 
utilizing  the  9400  Sylvania  computer.  The  necessary  computer  routine  to  accom¬ 
plish  this  evaluation  is  represented  by  the  flow  diagram  shown  in  Figure  24. 

Table  16  describes  the  underlying  frequency  distribution  of  the  various 

conriponents  in  the  circuits.  The  p  represents  the  nominal  value,  the  a  repre- 

2 

sents  the  standard  deviation,  (o  =  variance),  the  upper  and  lower  limits  repre¬ 
sent  the  truncation  points  of  the  distribution,  i.  e.  ,  the  tolerance  limits  for  each 
component. 

These  underlying  frequency  distributions  which  define  the  parameters  in  a 
statistical  manner  are  utilized  in  the  repeated  computation  of  the  transfer  func¬ 
tion.  The  results  of  1000  iterations  are  shown  in  Figures  25  through  28.  Having 
determined  acceptable  boundary  conditions  for  each  of  the  critical  performance 
criteria  required  for  acceptable  operation  of  the  circuit,  it  i»  possible  to  deter¬ 
mine  from  Figures  25  through  28  the  over-allreliability  of  the  circuit.  Moreover, 
the  sensitivity  of  the  individual  performance  criteria  as  a  function  of  the  varia¬ 
bility  of  the  individual  circuit  components  is  also  established.  This  process  can 
be  repeated  for  any  number  of  hours  of  simulated  use.  This  is  accomplished  by 
simulating  the  aging  of  component  parts  using  degradation  rates  as  shown  in 
Table  16.  These  rates  must  be  determined  by  analyzing  the  physics  of  failure  of 
various  elements  of  a  circuit.  [20] 

Referring  to  Table  16,  the  means  and  variances  (  p.  s  and  a  s)  for  the 
various  components  (at  t  =  1000  hours)  are  given.  With  this  information  that 
describes  the  underlying  frequency  distribution  of  the  parts,  it  is  possible  to 
determine  the  cumulative  distribution  functions  of  the  circuit  parameters  at 
t  =  1000  hours  by  repeating  the  technique  outlined  above.  The  cumulative  dis¬ 
tribution  functions  are  plotted  in  Figures  29  through  34.  As  a  check  on  the  ran¬ 
domness  of  the  number  generator,  two  of  the  circuit  input  cumulative  distribution 
functions  were  plotted.  At  t  =  0,  the  cumulative  distribution  function  of  Rj  was 
plotted  and  is  presented  in  Figure  35.  At  t  =  1000  hours,  the  cumulative  distri¬ 
bution  fxmction  of  was  plotted  and  is  shown  in  Figure  36.  The  general  shape 
of  these  curves  i.  e. ,  the  straight  line  indicating  uniform  distribution  and  the  S 
shape  curve  indicating  normal  distribution,  validate  the  correct  operation  of  the 
random  number  generator. 
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TABLE  16 


UNDERLYING  A  FREQUENCY  DISTRIBUTION  OF  PARAMETERS  USED  IN 
LOW-LEVEL  LIMITING  AMPLIFIER 


t'O  HRS  t-  1000  HRS 


MEAN 

P- 

STANDARD 

DEVIATION 

<r 

LOWER 

LIMIT 

UPPER 

LIMIT 

R1 

1,000 

29 

950 

1,050 

R2 

10,000 

290 

9,500 

10,500 

R3 

220 

6 

209 

231 

R4 

5,600 

162 

5,320 

5,880 

R5 

1,500 

43 

1,425 

1,575 

R6 

1,500 

43 

1,425 

1,575 

R7 

2,200 

64 

2,090 

2,310 

R8 

22,000 

640 

20,900 

23,100 

R9 

470 

14 

446 

493 

RIO 

330 

10 

313 

346 

Rll 

^100,000 

2890 

95,000 

105,000 

•pi 

20 

4 

10 

38 

•pz 

20 

4 

10 

38 

rbl 

10,000 

2890 

5,000 

15,000 

rb2 

10,000 

2890 

5,000 

15,000 

•Pi 

75 

26 

30 

120 

MEAN 

P 

STANDARD 

DEVIATION 

a 

LOWER 

LIMIT 

UPPER 

LIMIT 

1,000 

58 

900 

1,100 

10,000 

578 

9,000 

11,000 

220 

13 

198 

242 

5,600 

324 

5,040 

6,160 

1,500 

87 

1,350 

1,650 

1,500 

87 

1,350 

1,650 

2,200 

127 

1,980 

2,420 

22,000 

1270 

19,800 

24,200 

470 

27 

423 

517 

330 

19 

297 

363 

100,000 

5780 

90,000 

110,000 

27 

7 

10 

50 

27 

7 

10 

SO 

10,000 

2890 

5,000 

15,000 

10,000 

2890 

5,000 

15,000 

75 

36 

25 

150 

NOTL  ALL  DISTRIBUTIONS  UNIFORM  EXCEPT  &  fiZ  WHICH  ARE  TRUNCATED  NORMAL 
*TEXAS  INSTRUMENT  TRANSISTOR  RELIABILITY  DATA  3RD  gUARIIR  1961 
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Figure  25.  Cumulative  Distribution  Function 
of  V^2  (Zero  Time) 


Figure  26.  Cumulative  Distribution  Function 
of  (Zero  Time) 


Figure  27.  Cumulative  Distrflsution  Function  Fl^ire  28.  Cuimilotlve  Distribution  Function 
of  (Zero  Time)  of  (Zero  Time) 
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I  FIgur*  31 .  Cumulotlv*  Dlitrlbutfon  Function  of  V^2 

?  * 


f 


65 


FREQUENCY 


Figure  312.  Cumulative  Distribution  Function  of 
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Ftgur*  34.  Cumulativs  Distribution  Function  of  (t  <■  1000  Hrs.) 


11.84 
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FIgura  35.  Cumuiatlv*  Dlitrlbution  Function  of  R|  (Zero  Thno) 


69 


D 


5 


20 


40 


25  30  35 

( t  ■  1000  HRS  ) 


Figure  36.  Cumulative  Distribution  Function  of  ^ 


»  1000  Hrs.) 


As  an  example  of  the  usefulness  of  this  technique,  consider  the  critical 
performance  limit  of  ^E3  to  be  0.  6  volts.  Referring  to  Figure  37,  the  cumula¬ 
tive  distribution  function  of  at  t  -  0,  it  may  be  seen  that  a  0.  6  critical  per¬ 

formance  limit  produces  no  failures.  However,  making  use  of  the  Monte  Carlo 
technique,  referring  to  Figure  32  it  may  be  seen  that  at  t  =  1000  hours  that  4  per¬ 
cent  of  the  units  would  fail.  Thus,  the  example  above  briefly  highlights  the  use¬ 
fulness  of  this  technique  in  predicting  the  future  reliability  of  a  circuit  or  system. 

3.  2.  5  Transfer  Functions  of  the  Monostable  Multivibrator 

The  monostable  multivibrator  circuit  appearing  on  page  55  of  the  RADC 
Technical  Report  RADC-TR-59-243  has  been  analyzed  and  the  following  transfer 
functions  as  a  function  of  individual  circuit  parameters  as  well  as  the  supply  volt¬ 
ages  have  been  developed: 

(a)  Quiescent  current  and  voltage 

(b)  Output  amplitude 

(c)  Output  pulse  width 

The  schematic  of  the  monostable  multivibrator  is  shown  in  Figure  38  and 
the  performance  criteria  of  this  circuit  is  presented  in  Figure  39. 

The  monostable  multivibrator  has  two  states,  one  permanently  stable  and 
one  quasi-stable  state.  This  type  of  multivibrator  requires  a  triggering  signal  to 
change  from  the  stable  to  the  quasi- stable  state.  It  is  possible  for  the  multivibra¬ 
tor  to  remain  in  its  quasi- stable  state  for  a  long  period  of  time  in  comparison  to 
the  time  required  for  transition  between  states.  However,  no  external  signal  is 
required  to  reverse  this  transition,i.  e.  ,  eventually  the  multivibrator  will  return 
from  the  quasi- stable  state  to  its  stable  state  unaided.  [24] 

When  a  single  negative  input  trigger  is  fed  to  the  monostable  multivibrator 
represented  in  Figure  38,  an  output  gate  with  a  controlled  width  and  a  fixed  ampli¬ 
tude  is  delivered  to  its  output  terminal.  Transistors  and  Q2  form  a  regen¬ 
erative  feedback  pair  while  transistors  and  operate  as  emitter  followers 
providing  a  low  output  impedance  as  well  as  isolating  the  timing  and  trigger  func¬ 
tions  from  external  disturbances.  Transistor  Q2  is  normally  operating  in  the 
non-saturated  condition  while  is  cut  off.  When  a  trigger  is  fed  to  the  multi¬ 
vibrator,  Q2  is  cutoff  and  regenerative  action  drives  into  saturation.  During 
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FREQUENCY 


Figure  37.  Cumulative  Diitribution  Function  of  at  Time  Zero 
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TERMINAL 


Not*;  All  RMistors  +  5%  Toleronc*,  1/2  watt  Rating. 


Figure  38.  Schematic  of  Monoitable  Multivibrator 
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Figure  39.  Performance  Criteria  of  Monoftobie  Multivibrator 
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,  the  transition,  a  charge  is  stored  on  the  parallel  combination  of  capacitor  and 

the  external  timing  capacitor.  The  charge  on  these  capacitors  dissipates  in  a 
finite  time  period  corresponding  to  the  gate  length. 

The  basic  application,  that  of  a  gating  or  timing  circuit,  of  this  type  of 
multivibrator  results  from  the  fact  that  it  may  be  used  to  establish  a  fixed  time 
interval,  the  beginning  and  end  of  which  are  marked  by  an  abrupt  discontinuity  in 
the  voltage  waveform. 

Quiescent  Current  and  Voltage  for  Off  During  Stable  State 


The  DC  quiescent  currents  and  voltages  are  derived  when  the  monostable 
multivibrator  is  in  its  stable  state  with  no  input  signal  applied.  When  the  circuit 
is  in  this  condition,  Q^,  CR^,  and  CR^  can  be  removed  because  they  are  reverse 
biased.  Then,  Figure  40  is  the  DC  equivalent  circuit  for  the  state  in  which  Q^  is 
nonconducting  and  is  conducting.  Transistors  and  are  emitter 
followers  and,  therefore,  are  conducting  but  are  not  saturated  nor  cut  off.  A  0,  6  V 
drop  was  assumed  across  each  transistor  input  and  also  across  CR^,  CR^i  ^^4 
and  CRg  since  they  are  forward  biased.  The  cathode  of  zener  diode  CR^  was 
assumed  to  be  at  or  less  than  6.  2V  because  of  the  zener  diode  eqmvalent  circuit 
in  Figure  41.  A  complete  circuit  analysis  of  the  multivibrator  for  this  state  yields 
equations  (89)  through  (99). 


25  -  Vj.3  -  0.  6 


-(25  -  Vj.3)  ^  Vj.3  -  _ _ 

Ri  iq -  "  TT 


X3 


(89) 


^C3  -JT 


Ve3  "  25 


LI 


(90) 


Ve3-Ve4-1-2  .Ic4,  Ve4+1-2“Vi  25.Vj,4.0.6  ,  ^ 

iq - =-r  ^  - - ^  ^C2 


(91) 


2^-^E2-^-2  _  Vj,2+1.2  Ic2 

irrmrrr  “  irmrr,  ir  it; 


14  '15 


‘16  +  ‘'12 


‘8 


(92) 


I 

% 


5- 


75 


Ve2  -  ^1  .  Ve4  "  ' 

.2  - 

h  .  ^ 

(93) 

^8  ^7 

+  ^C2  ,  ^E2 

^cz  *  Tl 

(94) 

I  V 

I  +  _£i  = 

•f 

Ve4  -  ^2 

(95) 

P4  iqj 

^10 

Ve4  -  ^2  _  V2  +  25 

(96) 

''cz  *  ''e4  * 

(97) 

Ve4  "  ‘-2  - 

Due  to  CR6 

(98) 

''cz  ’  ''e2  * 

For 

Q2  Non-Saturated 

(99) 

These  equations  completely  describe  the  quiescent  currents  and  voltages  of  the 
monostable  multivibrator.  It  must  be  noted  that  in  order  for  the  transistor  to 
be  saturated  the  inequalities  (98)  and  (99)  must  be  satisfied.  Therefore,  it  is 
poss^ible  to  obtain  four  different  solutions  when  Equations  (89)  through  (97)  are 
solved.  The  following  procedure  will  indicate  which  of  the  solutions  is  correct. 

(a)  If  the  solution  is  such  that  =  5V  and  V^2  ^  ^^,2  +  0. 2  then 

equations  (89)  through  (97)  are  valid. 

(b)  If  the  solution  is  such  that  $  5V  and  V^2  ^  ^E2  ^ 

Q2  is  saturated  and  equations  (89)  through  (97)  must  be  revised  as 
follows:  In  equations  (92)  and  (94)  set  ^^B2  equation  (100) 

Vc2*Ve2+0-2  UOO) 

(c)  If  the  solution  is  such  that  >  5V  and  V^2  >  ^ 

CRg  is  in  zener  breakdown,  Q2  is  not  saturated  and  equations  (89) 
through  (97)  must  be  revised  as  follows: 
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In  equations  (89),  (91),  (93),  (95),  (96)  and  (97)  set  a  5V. 

(d)  If  the  solutionis  such  that  ^  5V  and  “  ^EZ  ^ 

is  in  zener  breakdown,  is  saturated  and  the  equations  must  be 

revised  as  in  (b)  and  (c). 

Equations  (89)  through  (97)  are  now  used  to  form  the  matrix  equation  (101). 
It  now  must  be  determined  which  of  the  four  conditions  listed  above  is  appropriate. 
Table  17  presents  the  computer  solution  of  the  matrix  equation  (101).  As  may  be 
seen  from  the  solution,  the  inequality  of  equation  (98)  is  not  satisfied.  Therefore, 
matrix  equation  (101)  is  not  appropriate  and  condition(c)  must  be  used. 

TABLE  17 

COMPUTER  SOLUTION  OF  MONOSTABLE  MULTIVIBRATOR  FOR  OFF 


< 

m 

''e3 

^C2 

''l 

''2 

'C2 

'C3 

'C4 

4,67 

WBSM 

9.85 

10.45 

4.55 

1.54 

.0055 

.0096 

.0055 

Table  18  presents  the  computer  solution  of  the  multivibrator  utilizing 
condition  (c)  and  the  actual  measurements  made  on  a  breadboard  circuit  in  the 
laboratory.  As  may  be  seen  from  Table  18,  all  equalities  are  satisfied  and  the 
close  comparison  of  calculated  and  measured  values  validate  the  accuracy  of  the 
transfer  function. 

TABLE  18 

COMPARISON  OF  COMPUTED  AND  MEASURED  VALUES  OF  MONOSTABLE 
MULTIVIBRATOR  FOR  Qj  OFF  AND  CONDITION  C 


''e2 

''e3 

''e4 

''C2 

''l 

^^2 

'C2 

'C3 

'C4 

COMPUTED 

3.55 

16.91 

5.0 

5.6 

3.04 

-2.16 

0.00417 

0.0067 

0.0013 

MEASURED 

3.9 

15.7 

4.5 

5.2 

3.2 

-2.37 

— - 

— 

— 
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Quiescent  Current  and  Voltage  for  Q2  Off 


When  the  monostable  multivibrator  receives  a  negative  trigger  of  sufficient 
amplitude  to  change  its  State,  the  equivalent  circuit  of  Figure  42  with  Q2  off 
applies  at  the  beginning  of  this  quasi-stable  interval.  In  this  state,  Q2  , 
and  CR2  &re  out  of  the  circuit  because  they  are  reverse  biased.  Transistors 
Qj  and  Q^  are  conducting  in  the  active  region.  The  drop  across  the  base 
emitter  junctions  of  Q^,  and  and  across  CRj,  CR^  ,  CRg  and  CR.^ 

was  assumed  to  be  0.  6V  because  they  are  forward  biased.  The  cathode  of  zener 
diode  CR/  is  at  or  below  6.  2V. 


•‘f  If  I'KtO' 


(110) 


^2  ■  ^E1  ■  _  ^C1 


V^3  +  1.2  §  6.2 


Vci  >  +  0.2 


(111) 

(112) 


Equations  (103)  through  (110)  define  the  bias  voltages  and  currents  assuming 
is  not  saturated  at  the  start  of  the  quasi- stable  period.  Inequalities  (111)  and 
(112)  must  also  be  satisfied  or  a  procedure  similar  to  that  previously  outlined 
must  be  followed. 

Equations  (103)  through  (110)  are  used  to  form  matrix  equation  (113). 

The  matrix  is  then  solved  by  a  computer  and  .Table  19  presents  this  solution.  As 
may  be  seen  from  Table  19  the  inequalities  of  equations  (111)  and  (112)  are  not 
satisfied;  therefore,  the  matrix  equation  (113)  is  not  appropriate  and  again  condi¬ 
tion  (c)  must  be  used. 


TABLE  19 

COMPUTER  SOLUTION  OF  MONOSTABLE  MULTIVIBRATOR  FOR  Q2  OFF 


''ei 

''e3 

''e4 

''l 

''2 

''3 

''ci 

'C3 

'C4 

7.099 

20.13 

23.09 

9.40 

10.43 

10.21 

.00833 

.01015 

,01147 

Table  20  presents  the  computer  solution  utilizing  condition  (c)  and  the  laboratory 

measureuients  made  on  a  breadboard  of  the  circuit.  Thus,  the  accuracy  of  the 
transfer  function  is  determined. 

Output  Amplitude 

The  output  amplitude  is  the  difference  between  the  outputs  (Vg.^  or  Vg,^) 
computed  from  Figure  40  and  from  Figure  42.  The  measured  amplitudes  were 
AVg.3  a  16V  and  =  I6V  for  =  1000  micro-farad  and  pulse  width  » 

30  microseconds. 
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TABLE  20 


COMPARISON  OF  COMPUTER  AND  MEASURED  VALUES  OF  MONOSTABLE 
MULTIVIBRATOR  FOR  Q2  OFF  AND  CONDITION  C 


'^El 

''e3 

''e4 

''l 

''2 

''3 

'ci 

'C3 

'C4 

COMPUTED 

6.8 

5.0 

17.0 

3.8 

7,1 

5,2 

0.008 

0.004 

0.008 

MEASURED 

6.5 

4.0 

15.0 

... 

... 

... 

... 

... 

Monostable  PuIbc  Width 

The  time  that  the  monostafcle  will  remain  in  its  quasi- stable  state  is  a 
function  of  the  circuit  parameters  and  supply  voltages.  Figure  43  shows  the 
complete  circuit  and  its  similaritv  to  a  simple  monostable  with  the  addition  of 
emitter  followers  and  clamp  diodes.  The  voltage  at  the  junction  of  and 

can  be  determined  by  the  same  techniques  used  on  similar  circuits.  [21]  The 
collector  voltage  of  Qj  and  the  bise  voltage  of  are  constrained  to  two 
voltage  levels  by  means  of  diodes  CR^  and  CR^.  These  voltages  have  been  com¬ 
puted  previously  and  they  appear  at  V^- ^  shifted  by  approximately  0.  6V  at  a 
relatively  low  impedance.  An  eqx  ivalent  circuit  is  shown  in  Figure  44  for 
and  CRg  conducting. 


The  voltage  across  Rj2  steady  state  is: 

R, 


12 


=  (Vp2 


T2 


12 


16 


(114) 


However,  during  the  quasi- stable  state,  Q2  and  CRg  will  be  cut  off  and  the 
equivalent  circuit  will  appear  as  shown  in  Figure  45.  Applying  Thevenin's 
Theorem,  the  circuit  in  Figure  45a  can  be  transformed  to  that  shown  in  Figure 
45b,  where; 


R 


(R’h  +  ^16>  ^12 

*^12  '■  ^14  ■*■^16 


(115) 


% 


84 


i 

> 

I 


c  »  c .  +  c 

4  X 


V  = 


*'12 

k,,  R',  .  +  R 


25 


12 


14  "  *'16] 


(116) 

(117) 


FIgurt  45.  Thaywifn  Equivalent  Circuit  for  Monoctabla  Pulw  Width 


The  voltage  at  the  junction  of  C  and  R  will  appear  as  shown  in  Figure  46. 


Pigwredd.  Output  Vultags  Wovofonii 


where:  j— 

E  .  V  -  ..2,.,j_^  - 

and  the  capacitor  C  discharges  toward  V  so  that  the  voltage  at  the  junction  of 
R  and  C  during  the  quasi-stable  state  is: 


(Ve3  Off) 


E3 


■ 


(118) 


V  -  E  exp  (^) 


(119) 


12 

The  quasi-stable  state  will  terminate  when  e  =  (^E2  *  ^)  Ti - TTR — 

Therefore,  the  width  of  the  pulse  t^^  can  be  obtained  by  solving  the  following 
equation  for  t^. 


<''e2*  =  V.E  ,xp 


(^) 


^  **P  ^  ■  <^E2  f-  ^)  R~TTl^^ 


exp 


12 


V  -  (V-,  +  1.2)  w- . V-v.  - 


(120) 


Solving  for  t 


therefore. 


V 

«*P  ife  = 


_ E _  _ 

'Ve2*‘-^’(r,2+"S,4) 


t  *  RC  In 
P 


E 


V-(Vj.2+l 


(121) 


•T 


where: 


E  =  V-[(V£2  -r  1.  2)  -  On  -  Of£)l 


E3 


'E3 


and 


V  = 


(25) 


For  nominal  values  of  circuit  parameters: 


(122) 

(123) 


V  =  19.5V 

E  =  19.  5  -  [(3.  56  +  1.  2)  -  (16.  91  -  5.  0)1  =  25.  65V 

-  °  <55.^K)Clnl.74 

19.  5  -  4.  76 

tp  =  (31.  0  X  10^)  C  (124) 

Therefore,  it  may  be  seen  from  equation  (124)  that  the  pulse  width  is  directly 
proportional  to  the  value  of  the  external  capacitor  and  .  Table  21  pre¬ 
sents  a  comparison  of  the  computed  and  measured  values  of  the  pulse  width.  It 
may  be  seen  by  referring  to  the  table  that  the  values  check  closely,  thus  indicating 
the  accuracy  of  the  transfer  function! 


TABLE  21 

COMPARISON  OF  COMPUTED  AND  MEASURED  VALUES  OF  PULSE  WIDTH 


c 

X 

c  +C.+C^ 

X  4  stray 

Computed  t 

P 

MocMurod  t 

P 

PRF 

0 

50  pf 

1.56  mlcTOMc. 

2  mIcroMC. 

50  KC 

160  pf 

210  pf 

6.5  mleraNC. 

7  mieroMc. 

50  KC 

100  pf 

1050  pf 

32.5  mlcratoc. 

30  mlerasoc. 

5ICC 

.l/if 

.Ipf 

3.1  mllltioe. 

2.6  mlllhoc. 

SOcpi 

.Spf 

.jpf 

15.5  mllllMc. 

12.5  milllMe. 

aospi 
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3.  2.  6  Transfer  Functions  of  the  Bistable  Multivibrator 

The  bistable  multivibrator  circuit  shown  in  Figure  47  has  been  studied  and 
the  following  transfer  functions  have  been  developed; 

a.  Quiescent  current  and  voltage  as  a  function  of  individual  circuit 
parameters  as  well  as  the  DC  supply  voltage. 

b.  Output  amplitude  as  a  function  of  individual  circuit  parameters  as  well 
as  the  DC  supply  voltage. 

c.  Minimum  trigger  amplitude  as  a  function  of  individual  circuit 
parameter. 

Table  22  lists  the  performance  criteria  for  this  circuit. 


♦MV 


Note:  This  circuit  is  contained  in  RADC  TR  59-243,  dated  December  15, 
1959,  Titled  Reliable  Preferred  Solid  State  Functional  Divisions. 

Figure  47,  Bistable  MultIvflNvtor  SdMmatle 
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TABLE  22 


PERFORMANCE  aiTERIA  OF  THE  BISTABLE  MULTIVIBRATOR 


Roquiromonts 

Maximum 

Dosign 

■hmSiSmbo 

■HolnOOBlLoB 

Triggor  ^‘lip-Flop)  Input 

- 

- 

«p 

AC  Couplod 

m 

- 

- 

Triggor  Arnplltudo 

-8  V 

-5  V 

-4v 

Triggor  Rlsotimo 

1  |iSOC 

-  - 

0.1  psoc 

Triggor  Froquoncy 

250  kc 

- 

- 

Input  Impodonco 

- 

m 

9,ooon 

Sot-Rosot  (Bistoblo)  Input 

m 

m 

- 

DC  Couplod 

m 

- 

- 

Veitogo  Lovol 

+6  V 

+2  V 

0  V 

Rlsotimo 

- 

- 

0.1  psoc 

Pulso  Froquoncy 

250  kc 

- 

- 

input  Impodonco 

- 

- 

9,OOOQ 

Output  Arnplltudo,  Pook-to-Pook 

18  V 

15  V 

10  V 

Output  Polarity,  Both  DC  Couplod 

m 

- 

- 

Output  Rlsotimo 

1  psoc 

- 

0.1  psoc 

Output  Docoy  Timo 

1  psoc 

■0 

0.1 psoc 

Output  Loading 

- 

- 

2,500Q 

DC  Supply  Voltogos 

- 

+25  V 

- 

- 

-25  V 

- 

Oporating  Tomporoturo 

8^C 

- 

-55®C 

NOTEt 

This  Is  oontaliMd  In  RAOC  TR  59-243,  dotad  Dtcmihtr  15,  1959, 
TItInd,  Rnllobb  PrsfwTsd  Solid  Skdo  FwneHenai  Dtvblens,  Cen- 

tioct  AF  85' W1WV - ^ - * - 


Quiescent  Currents  and  Voltages 


Since  this  circuit  is  symmetrical,  the  DC  bias  conditions  will  be  computed 
for  the  bistable  multivibrator  circuit  in  one  state  (Q^  non-conducting,  Qj  conduct¬ 
ing).  Thus,  the  bias  will  be  known  for  the  opposite  stable  state. 


Ffgur*  48,  DC  Equivalent  Circuit  of  llstablo  with  off 


Figure  48  is  the  DC  equivalent  circuit  for  the  state  in  which  is  non¬ 
conducting  and  is  conducting.  .  CR^  ,  CR2 ,  CR^,  and  CR^  are  removed 
frmn  the  circuit  because  they  are  reverse  biased.  Figure  9  is  the  DC  equivalent 
circuit  for  ^a  conducting,  non-saturated  transistor.  Diodes  CR^,  CR^,  CR^  and 
CRg  were  assumed  to  have  a  0.  6  volt  drop  between  anode  to  cathode  because  they 
are  forward  biased.  Figure  10  represents  the  DC  equivalent  circuit  of  a  transis¬ 
tor  in  saturation.  The  equivalent  circuits  presented  in  Figures  9  and  10  were 
utilized  in  solving  the  transfer  functions  for  the  condition  of  non-saturated 
and  saturated  respectively. 

An  analysis  of  the  circuit  yields  equations  (124)  through  (132).  These 
equations  are  then  simplified: 


Cl  _ 

“r  ■ 


(128) 
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mmm 


(131) 


I 


B3 


(132) 


The  matrix  equation  (133)  is  then  solved  using  a  computer  and  the  results  are 
presented  in  Table  23. 


TABLE  23 

COMPUTER  SOLUTION  FOR 
BISTABLE  MULTIVIBRATOR 
WITH  Q4  OFF,  Q,  SATURATED, 
RlI  *  ^^2  “  2500, p  «  25 


V, 

O.fOSf 

''2 

-2.04 

''ll 

4.3MI 

''12 

u.o 

''0 

4.7M0 

'ci 

0,0012 

'€2 

0.00022 

O.OOSItO 

■n 

-0.000m 

NOTE:  Waveforms  at  th*  minimum  trl|f«r  Uv«l:  dottud 

wavnform*  occur  when  Innufficicnt  input  it  proaont. 


Figure  49.  Bistable  Equivalent  Circuit  at  the 
Leading  Edge  of  the  Trigger  Pulse 


The  solution  for  in  saturation  appears  satisfactory.  However,  a 
negative  answer  for  base  current  1^^  is  not  physically  possible.  The  negative 
answer  indicates  that  the  circuit  does  not  allow  sufficient  base  current  to 
saturate  as  it  would  if  the  solution  for  base  current  were  positive.  This 
means  that  the  assumed  condition  (Q^  saturated)  cannot  be  physically  realized. 

Therefore,  the  analysis  must  be  repeated  using  Figure  9  to  represent 
which  is  apparently  conducting  but  is  not  saturated.  This  change  results  la  the 
following  changes  in  the  equations  (124)  to  (132). 


(  This  page  intentionally  left  blank  ) 


in  equations  (124)  and  (132);  and  equation 


is  substituted  for 
(131)  is  eliminated. 

The  revised  matrix  equation  (134)  is  then  solved  and  the  results  together 
with  the  measured  results  are  presented  in  Table  24. 

TABLE  24 

BISTABLE  MULTIVIBRATOR  OFF  Qg  NONSATURATED 
p  =  25,  Rl,  =  Rl2=  2500 


''i 

''2 

''e2 

''e3 

'ci 

'C2 

Canputadt 

I.0S8  V 

-t.244V 

4.1W  V 

17.40  V 

4.014  V 

3.907  MA 

9.31  MA 

3.84  MA 

1,07  V 

-1.33  V 

4.1  V 

17.0  V 

4,2  V 

•M 

Coi«puM:  •  Vj,  +  0.4  -  S.80  V 

MtoMndi  V^*4.8V 


Table  25  shows  the  solutions  for  the  8  variables  for  transistor  p's  of 

10,  20,  .  .  .  60.  This  data  shows  the  changes  in  bias  due  to  changes  in  transis> 

tor  current  gain  and  indicates  that  the  output  amplitude  (V^2  ~  increased  by 

11.  4  percent  for  a  change  in  p  from  10  to  60. 


TABLE  25 

BISTABLE  MULTIVIBRATOR  OFF  Qg  NONSATURATED 
R^^  *  R|^2  *  2500,  “10,  20,  •  •  • ,  60 


''l 

'^2 

''ei 

n 

''d 

•ci 

'C2 

•1.010 

4.904 

17.07 

3.92B 

0.008SI9 

O.OOS719 

o.oinai7 

-to 

■■DHI/fli 

1 

i 

■KULW 

OB 

OsindvpQ 

_ 

■■DftZ'Sl 

■EKM 

HULliaJHI 

■C2  J 

ilFTC-H 

HKAiuJHi 

MCLJI 

■  II  mm 

—lEIMi 

HCKI3B 

■LULJBI 

HKJLjUHI 

■JLitiziM 

MC9LJI 

ItK'HMi 

■xim 

■ULliJBl 

■I'Kv.-iilJi 

■KXJi 

FIguni  50.  Equivalant  Circuit  for 
tho  Input  of  Qj  Prior  to  Tom  Off 


®c«i  •  "Sw  •  *«»  pf)  •  in  pf 
C|  >  109  pi 


Figure  51.  Equlvolof^  Ckeult  for  and 
tho  Input  ClKutt  Prior  to  Turn  Off  wMi 

WCUvtVM  TrWwiOnm 
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Table  26  is  the  computed  solution  of  the  bias  equations  assuming  no 
output  loading,  i.  e. ,  ~  ^  major  difference  between  this  case 

and  that  of  Table  24  is  that  the  collector  currents  of  and  are  reduced 
by  300  percent.  The  other  bias  conditions  changed  less  than  10  percent. 


TABLE  26 

COMPUTER  SOLUTION  FOR  BISTABLE  MULTIVIBRATOR  WITH  OFF  Qg  NONSATURATED 

p  s  p  s  eo 
•'Ll  ''L2 


''l 

''2 

Vei 

''e2 

''e3 

'ci 

'C2 

b 

Cempuladi 

1.08  V 

-1.15 

4.32 

17,7 

4.14 

0.00152 

0.00264 

0.00398 

MaoMiradt 

1.07  V 

-1.33 

6.1 

17.0 

4.2 

— 

— 

— 

Solution  for  Output  Amplitude 

The  output  amplitude  of  the  bistable  multivibrator  providing  sufficient 
trigger  voltage  is  present  to  change  its  state,  equals  the  difference  in  quiescent 
values  of  the  emitter  voltages  and  Q^,  i.  e. ,  the  output  amplitude  s 
(Vj,2  -  value  computed  from  Figure  12  using  nominal  values  for  all 

components  is  17.7  -  6.  3  =  11.  4V.  The  measured  value  obtained  irom  the  bread¬ 
board  circuit  was  IIV.  The  measured  and  calculated  values  compare  favorably, 
therefore,  confirming  the  accuracy  of  the  solution. 


Minimum  Trigger  Amplitude  , 

To  find  the  minimum  required  trigger  amplitude,  the  equivalent  circuits 
of  Figures  49,  50,  and  51  must  be  used.  From  Figure  49,  it  can  be  seen  that 
Vq  must  be  changed  negatively  by  at  least  0.  6V  to  nullify  its  0.  6V  "on"  bias. 

Then  will  rise  enough  to  propagate  through  and  turn  oh.  Figure 

50  shows  the  calculations  for  obtaining  the  input  capacitance  and 

resistance  of  Qg.  Figure  51  is  used  to  compute  the  required  to 

turn  Qg  off.  At  the  leading  edge  of  the  pulse,  from  the  circuit  Figure  51: 


HA 


/ 


* 

/ 


''bE3  ” 


S^BE3'^=.I 

.  "I: - 1 - ■CR4 

^5  •l=BE3+C4> 


V™,  .  V.  /, 


‘”'=5^  W« 


r)-' 


CR4 
(135) 


/Cc  +  C  e.  +  C  \ 

The  required  =  ^^BE3  ^  '^CR4^  ( - C - )  terms  of  the 

circuit  parameters  and  specified  transistor  parameters, 


V.  = 
in 


*^BE3'‘'^CR4^ 


1 


S  ^  2it(p  +  1)(0.26/Ij.)f«  ■*■  Cl“^ 


^2  P^oB2 


Substituting  nominal  values  into  the  above  equation  yields 
V.^  =  [150+  200+  63.5]  =  3. 3V 


(136) 


(137) 


The  measured  minimum  required  trigger  amplitude  was  3.  5V  peak  to  peak 
at  the  "trigger"  input.  This  input  becomes  differentiated  by  Cg  and  R^^  so  that 
a  negative  spike  of  1.  5V  amplitude  occurs  at  the  junction  of  Cg  and  R^^  when 
the  input  is  going  negative.  This  appears  as  a  negative  step  of  0.  8V  amplitude  at 
This  step  instantaneously  raises  the  collector  of  by  3V  and  is  propaga¬ 
ted  through  the  emitter  follower  to  the  base  of  Q^.  It  appears  as  a  positive 
step  of  2.  IV  at  thus  raising  the  base  of  to  +0.  8V  thereby,  causing 
to  conduct.  The  negative  voltage  V^^  is  fed  back  through  Q2  to  the  base  of 
turning  off. 

Thus,  the  procedure  described  above  details  the  stepslused  in  the  devel¬ 
opment  of  three  transfer  functions  for  the  bistable  multivibrator  circuit. 

3.  2.  7  Transfer  Functions  of  the  Pulse  Adder 

The  pulse  adder  analysed  is  basically  a  form  of  a  wide  band  video  amplifier. 
The  basic  amplifier  schematic  is  shown  in  Figure  52.  This  sunpliRer  can  be  used 
to  perform  the  function  on  an  OR  circuit,  which  is  basically  a  mixer  or  buffer  that 
permits  a  number  of  pulses  to  be  connected  to  a  cmunon  load.  It  also  minimises 
the  interaction  of  the  pulse  sources  on  each  other.  If  a  positive  pulse  with  a 
maximum  amplitude  of  5  volts  is  applied  to  either  of  the  ii^ut  terminals,  a  palse 


Figure  52.  Pulse  Adder  Schenratlc 

of  similar  amplitude  and  form  will  be  available  at  the  output.  However,  if  a  pulse 
is  applied  at  both  terminals  simultaneously,  the  output  consists  of  the  sum  of  the 
two.  If  the  pulses  are  of  sufficient  amplitude  to  saturate  the  stage,  then  only  the 
wider  of  the  two  is  transmitted,  i.  e. ,  a  pulse  of  sufficient  width  and  amplitude 
inhibits  and  prevents  the  transmission  of  a  smaller  or  narrower  pulse  through  the 
circuit.  The  performance  criteria  of  this  circuit  is  presented  in  Table  27. 


TABLE  27 

PERFORMANCE  aiTERiA  OF  PULSE  ADDER 


mamm 

Moximum 

Design 

Cmiter 

Minimum 

Input!  (Two  Unboloneod 
to  Ground) 

Input  Impodanco 

— 

war 

8000 

Inputs  Potitivo  Paok 

5v 

3.5  V 

— 

Source  Impodanco 

750 

— 

Voltogo  Gain,  Either  input 

1 

0.75 

0.6 

Output  Impodanco 

~ 

750 

~ 

Bandwidth 

4  me 

mm 

lOcps 

RIm  Time 

O.iS/jsec 

0.1  piec 

-- 

Output  (One  Unbalanced) 

mm 

MO! 

— 

Output  Level 

5v 

3.5  V 

»  _ 

Output  Polarity  Poeltive 

Repetition  Rote 

50(X)  ppt 

-- 

20ppe 

Pulie  Duty  Cycle 

0.20 

— 

— 

Additive-Factor 

0.25 

load  Impedance 

— 

750 

DC  Supply  Voltage 

mm 

+12  V 

— 

Operating  Temperature 

85*C 
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Figure  53.  DC  Equivalent  Circuit  of  PuIm  Adder 


The  equatione  are  then  eimplified  and  are  used  to  form  matrix  equation  ( 144). 
Then,  aseuming  all  components  to  be  at  their  nominal  value,  p  ^  *  ^2  * 

>75  ohms,  the  matrix  equation  is  then  solved  using  a  computer.  Also,  a 
breadboard  circuit  of  this  pulse  adder  was  constructed  and  tested  in  the  laboratory. 
Table  28  presents  a  comparison  of  the  measured  and  calculated  values  of  the 
circuit  parameters.  As  nuny  be  seen  from  Table  28,  the  calculated  and  measured 
values  c<mpare  quite  favorably,  thus  validating  the  accuracy  of  the  transfer 
functiem. 
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TABLE  28 


COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF  THE  PULSE  ADDER 


N 

''bi 

Vb2 

''e2 

''c 

•ci 

•C2 

Calculated 

• 

0.98 

0.98 

0.38 

11.57 

Meanired 

0.93 

0.94 

0.30 

m 

— _ 

Maximum  Output  Amplitude 

The  maximum  peak  to  peak  output  amplitude  is  the  difference  between  the 
output  (Vg^)with  and/or  Q2  saturated,  and  the  output  (V^2)  u>'*^**‘ 
conditions.  The  voltage  y^2  l>een  previously  calculated  and  found  to  be  0.  38 
volts.  In  order  to  find  the  voltage  with  Q|  and/or  Q2  saturated,  the 

equivalent  circmt  of  Figure  54  must  be  used.  The  analysis  of  this  circuit  yields 
equations  (145)  through  (147|.  These  equations,  when  rearrange^  form  the  matrix 
equation  (148). 


12 


''E 

tt: 


0.  6 


■Vj,  -  0.  6 


+  I, 


(145) 


Vc  =  Vj.  +  0.2 


(146) 


12  -  V 


TT 


3 


(147) 
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Figur*  54.  Equlvalant  Circuit  of  PuIm  Addor  with  and/or  Q2  Soturotod 


r  1  ^  ^ 

-B —  +  XT"  *  TT" 
^2  *^4 

0 

-1 

’11.4  0.6" 

■pq"  •  iq* 

-1 

1 

0 

3 

0.2 

0 

1 

1 

12 

[_  ^  J 

(148) 


This  matrix  is  solved  using  a  computer,  and  the  results, along  with  those  made  on  a 
breadboard  circuit  are  presented  in  Table  29. 


TABLE  29 

COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF  THE  PULSE  AODBt 

WITH  Q,  an4/or  Q2  SATURATED 


''61 

''c 

*C 

CoieulotMli 

5.38 

5.57 

0.086 

Mooiuradt 

4.8 

5.0 

— — 
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Since  the  maximum  output  amplitude  ia  defined  ae  ^E1  ■  ^E2'  it  ia  poa  Bible  to 
determine  the  deaired  reaulta: 

Maximum  Amplitude  •  ■  5.  38  -  0.  38  «  5  volts 

The  above  analysia  allows  the  transfer  function  relating  input  and  output  voltages 
to  be  graphed  as  shown  in  Figure  55. 


3.  2.  7.  1  Determining  Incremental  AC  Gain  Utilizing  the  Signal  Flow  Graph 
Technique 

In  order  to  determine  the  incremental  AC  gain,  the  AC  equivalent  circiiit 
of  the  pulse  adder  (Figure  56)  must  be  constructed.  An  analysis  of  this  circuit 
yields  equations  (149)  through  (151).  The  signal  flow  graph  is  then  constructed 
from  the  equivalent  circviit  and  the  equations  and  is  presented  in  Figure  57. 
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+  I 
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^cl  ^  ^ 

^  ^’bl  ’^bl  ’^bl 


e. ,  -  0.  6  e 
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^c2  ^  ^  _  2°. 
^  ’^bZ  'bZ 


0.6  ^  *iZ  ~  ^ 

'bZ  ’^bZ 


(151) 


e 

It  is  now  possible  to  determine  the  incremental  voltage  gain  from  the 


flow  graph  and  in  equation  (152) 


in 
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lei 


Figur*  57.  Flow  Gioph  for  Pulso  Addor  Voltogo  Goto 


3.  2.  8  Transfer  Functions  of  the  Distribution  Amplifier 


The  distribution  amplifier  (emitter  follower)  functions  as  a  general  purpose 
amplifier  in  a  similar  manner  to  a  chathode  follower.  It  is  designed  to  accept 
short  duty  cycle  pulses  and  to  provide  power  gain  with  a  voltage  gain  of  mity*  The 
input  impedance  of  this  amplifier  is  normally  higher  than  the  load  resistance  of 
the  source.  Thus,  the  parallel  connection  of  several  amplifiers  to  a  cofamoa 
source  can  be  effected  without  loading  the  inptd.  Isolaticm  is  provided  hetw^as  in¬ 
put  and  output  signals. 
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The  schematic  of  the  amplifier  ie  shown  in  Figure  58  and  its  performance 
criteria  is  shown  in  Table  30.  The  ii^ut  impedance  of  this  transistor  amplifier  is 
frequency  sensitive  and  appears  reactive  above  the  p  cutoff,  frequency.  Compen> 
sating  for  this  is  the  network  LjR^R2  shown  in  Figure  58.  R^  and  R2  provide 
circuit  stability  by  determining  the  quiescent  operating  points  of  the  transistor. 
Diode  CRj  acts  as  a  direct  current  restorer  permitting  operation  with  duty  cycles 
up  to  50  percent.  Optimum  performance  is  achieved  by  restricting  duty  cycles 
to  20  percent. 


FlgiirsSB.  SdMmcdte  of  Hw  DMrllMitlon 
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TABLE  30 

PERFORMANCE  CRITERU  OF  DISTRIBUTION  AMPLIFIER 


Raqulramanh 

Maximum 

Design 

Center 

Minimum 

Input  (On«,  Unbalanced  to  Ground) 

— 

... 

Input  Impodonca 

— 

— 

soon 

Input  Lovol,  Potitivo  Poole 

5 

3.5  V 

— 

Source  Impedance 

. — 

750 

— 

Voltage  Gain 

1 

0.75 

0.6 

Bandwidth  (3db  Down  Frequencies) 

4  me 

— 

10  cps 

Rise  Time 

0.15  ^sec 

0.1  ^sec 

— 

Output  (One,  Unbalanced) 

— 

— 

— 

Output  Polarity  Positive 

Load  Impedance 

— 

750 

500 

DC  Supply  Voltage 

— 

+12  V 

— 

Operating  Temperature 

85®C 

— 

-55®C 

A  nodal  analysic  of  this  eqtdvalent  circuit  results  in  equations  (153)  and 


(155). 


12  - 


Vj,.0.6 


Vg  +  0.  6 


I 

c 


1 


c 


(153) 


(154) 


(155) 


These  equations  when  rearranged  form  the  matrix  equation  (156). 


'  1  ..  1 

1  ’ 

’11.4  0.6' 

ir/iq 

0 

F 

0 

1 

1 

V 

c 

m 

12 

^2 

1 

0 

I 

c 

0 

(156) 


The  matrix  shown  in  equation  (156)  was  programmed  on  the  Burroughs  205 
computer  and  the  solution  is  shown  in  Table  31. 


TABLE  31 

i  COMPUTSI  SOLUTION  OF  DISTRIBUTION  AMFLFIBI 


''t 

.453 

11.9 

.OOM 

lot 


The  input  voltage  required  to  saturate  the  transistor  can  be  obtained 
by  analyzing  the  equivalent  circuit  shown  in  Figure  60. 


Figure  60.  Equivalent  Circuit  of  Distribution  Amplifier  with  Saturated 

A  circuit  analysis  results  in  equations  ( 157)  to  ( l6l). 
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Vb  -  Vg  +  0.6 


(160) 


Vc  -  Vj,  +  0.2 


(161) 


These  equations  when  rearranged  form  the  matrix  equation  (162), 


-1 

0 

1 

0 

o" 

0.6 

0 

1.1,1 

0 

1 

0 

e  * 

A 

+  *2 

-1 

0 

1 

0 

0 

s 

0.2 

1 

0 

0 

1 

1 

0 

0 

0 

1 

^3 

0 

1 

• 

- r 

(162) 


(1.5,  3.5,  5.0,  7.0,  8.0,  9.0,  10. 0,  10.2,  10.4,  10.6,  10.8, 

11,  12) 


To  determine  the  minimum  input  voltage  required  to  saturate  the 
computer  was  programmed  to  search  for  the  value  of  e^  which  results  in  a 
negative  base  current  which  would  indicate  that  the  amplifier  is  in  saturation. 

The  results  of  this  computer  program  is  shown  in  Table  32.  This  shows  a  transi¬ 
tion  of  positive  to  negative  base  current  occurring  at  an  input  voltage  of  10.  677 
volts. 


3.  2.  8.  1  Small  Signal  Gain  of  Distribution  Amplifier 

To  find  the  small  signal  gain  the  equivalent  circuit  shown  inIFigure  61  is 
utilized. 

A  circuit  analysis  is  performed  and  equations  (163)  to  (167)  are  formula^rf. 


(161) 


^  ib  •  0 


(164) 
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TABLE  32 


DISTRIBUTION  AMPLIFIER  WITH  SATURATED 


''e 

''b 

'^c 

'b 

'c 

•fl 

5.02 

5.22 

0.0678 

1.5 

6.04 

6.24 

0.0576 

3.5 

6.81 

7.41 

7.01 

-0.0354 

0.0499 

5.0 

7.84 

8.44 

8.04 

-0.0229 

0.0396 

7.0 

8.36 

8.96 

8.56 

-0.0167 

0.0344 

8.0 

8.87 

9.47 

9.07 

-0.0104 

0.0293 

9.38 

9.98 

9.58 

-0.00421 

0.0242 

10.0 

9.49 

10.1 

9.69 

-0.00297 

0.0231 

10.2 

9.59 

mSM 

9.79 

-0.00172 

0.0221 

10.4 

9.69 

■M 

9.89 

-0.00048 

0.0211 

10.6 

9.72 

10.3 

9.92 

-0.00016 

0.020B 

10.65 

9.72 

10.3 

9.92 

-0.00010 

10.66 

9.72 

10.3 

9.92 

-0.00010 

0.0197 

10.67 

9.73 

10.3 

9.93 

-0.00002 

0.0207 

10.673 

9.73 

10.3 

9.93 

-0.00001 

0.0207 

10.674 

9.73 

10.3 

9.93 

0 

0.0207 

10.675 

9.73 

10.3 

9.93 

0 

0.0207 

10.676 

9.73 

10.3 

9.93 

0 

0.0207 

10.677 

9.73 

10.3 

9.93 

0.00001 

0.0207 

10.678 

9.73 

10.3 

9.93 

0.00002 

0.0207 

10.679 

9.73 

10.3 

9.93 

0.00002 

0.0207 

10.68 

9.74 

10.3 

9.94 

0.00006 

0.0206 

10.69 

9.74 

10.3 

9.94 

0.00015 

10.7 

9.77 

10.4 

9.97 

10.75 

9.79 

10.4 

9.99 

0.00077 
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10.5 

10.1 

0.0190 

11.0 

0.0199 


Using  equation  (163)  to  solve  for 


„  _  in  A 

-  R  rf  R." 

A  in 


R.R. 

A  in' 


(168) 


Using  equations  (165),  (166),  (167),  and  (168)  the  signal  flow  graph  shown  in  Figure 
62  can  be  obtained  and  the  gain  determined,  [22] 


Figure  62.  Signal  Flow  Graph  of  DhtrAMition  Amplifier 


Ra*^4(1  +  P) 


(169) 
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Where 


«  2200  (22.  000) 
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% 

*  50 
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4 

in 
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o 
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A  conventional  analysis  using  nodal  equations  gives  the  same  results.  Equations 
(172),  (172)  and  (174)  are  the  nodal  equations  for  the  equivalent  circuit  shown  in 
Figure  61.  This  is  shown  in  matrix  form  in  equation  (175). 
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e. 
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1 

0 
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o 
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Equation  (176)  shows  the  matrix  with  nominal  circuit  parameters,  and 
Table  33  is  the  solution  to  the  matrix  obtained  by  using  the  Burroughs  205 
computer. 


0.  013833  .  .  . 

1 

0 

0 

0. 0133 

0 

51 

-0. 002127 

^b 

= 

0 

-0.  02 

1 

0.  02 

e 

_  0  _ 

0 

(176) 


TABLE  33 

COMPUTER  SOLUTION  OF  DISTRIBUTION  AMPLIFIER  SMALL  SIGNAL  GAIN 


''b 

•b 

a 

0 

.96147V 

40 

.95947V 

3.  3  TRANSFER  FUNCTIONS  OF  MICROMINIATURE  CIRCUITS 

This  section  presents  the  transfer  functions  of  the  microminiature  circuits 
developed  by  the  Sylvania  Microelectronics  Laboratory. 

The  circuits  presented  in  this  section  of  the  report  were,  in  some  in¬ 
stances,  analyzed  not  simply  as  individual  circuits  with  their  supply  voltages> 
but  also  with  respect  to  their  use  in  the  system  represented  in  Figure  6.  An 
example  of  this  is  the  steering  circuit  which  was  analyzed  with  consideration  for 
the  loading  of  the  previous  stage.  Also,  the  active  circuits  (the  Flip-Flop  and 
NOR  circuit)  were  analyzed  with  their  respective  transistors  operating  in  both 
the  saturated  and  non- saturated  states.  The  accuracy  of  the  transfer  functions 
was  verified  by  measuring  the  voltages  on  actual  circuits  in  the  laboratory. 

Since  the  flip-flop  is  of  a  symmetrical  nature the  transformaticm  matrix 
solved  for  Q2  conducting  and  off  also  represents  the  condition  when  is 
conducting  and  Q2  is  off. 


« 


3.  3.  1  Transfer  Function  of  the  5  me  Flip-Flop 


The  flip-flop  or  bistable  multivibrator  circuit  utilized  for  analysis  is 
shown  in  Figure  63. 


4.SV  12V 


Figwe  63.  Five  Me.  Flip-Flop  Circuit 


This  general  purpose  type  circuit  is  specifically  designed  to  operate  in 
the  saturated  mode  at  a  maximum  repetition  rate  of  5  me.  This  type  of  circuit 
finds  extensive  application  in  pulse  circuitry  being  used  not  only  for  the  generation 
of  square  waves  from  pulses  but  also  for  the  performance  of  certain  digital 
operations  such  as  counting.  The  actual  operation  of  the  circuit  is  similar  to  the 
one  previously  discussed  in  Section  3.  2.  6. 

Quiescent  Currents  and  Voltages 

When  the  condition  exists  that  is  off  and  Q2  is  in  either  the  satuvitie4 
or  the  non- saturated  state,  Q I  and  CR2  and  the  associated  cireititry  art  re¬ 
moved  from  the  circuit.  Figure  64  is  the  equivalent  circuit  used  for  the  amdlyeitt 
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TABLE  34 

PERFORMANCE  CRITERIA  OF  MICROMINIATURE  FLIP  FLOP 


Moxlffluin  operating  frequency  *  5  me 

INPUT  REQUIREMENTSt 

Input  puiM  derived  from  pulra  level  gate 
Mfficimum  "Fen  in"  (puira  level  gates)  -  3  per  input 

OUTPUT  REQUIREMENTS: 

R|^  mtn  (max  bad  to  ground)  ■  6K  (0.8  ma) 

Rj  mln  (max  load  to  V^^)  «  2K  (6  mo) 

Maximum  DC  "Fan  out"  "  2  logic  inputs  or 

2  le^l  inputs  (j>ulse  level  gate) 
Maximum  AC  "Fan  out"  *  4  pulse  inputs  (pulse  level  gate) 
Maximum  capacitive  bad  »  50  pf 

Typical  output  clxiracteristics  for  operation  with  pulse 
level  gate  output  at  2j^C  unloaded 

turn  on  time  *  30  ns 
turn  off  time  ■  50  ns 
turn  on  delay  *  25  ns 
turn  off  delc^  »  20  ns 

Logic  levels:  "false"  *  -*-5  volts 
"true"  ■  +0.3  volts 

Supply  Voltoges 

V  +12  volts +5  percent 
cc  • 

''bb-<  volts  +5  percent 
> +4.5  volts +5  percent 
Vgg  ■  ground 

Storage  temperature  -  to  +  12a^C 
Operating  temperoture  ■  to  +  II^C 
Fewer  dtisipatien  ■  90  mw 
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Figura  64,  DC  Equivaivnt  Circuit  of  Flip*f  lop  with  off  and  Q2  Conducting 

The  equivalent  circuits  used  for  the  transistors  are  presented  in  Figures 
9  and  10.  Utilizing  the  above  set  of  equivalent  circuits,  the  analysis  yields 
equations  (177)  through  (18Z);  these  equations  are  then  rearranged  and  foma 
matrix  equation  (183). 
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This  matrix  equation  is  then  solved  using  a  computer.  The  results  of  this  com¬ 
puter  solution  along  \vith  the  results  received  from  laboratory  measurements  on 
this  circuit  are  presented  in  Table  35. 


TABLE  35 

COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF  FLIP-(=LOP 
WITH  OFF  AND  Qj  SATURATED 


CALCULATED 

MEASURED 

''cz 

^B2 

''ci 

im 

'c 

in 

0.78 

0.85 

5.1 

5.0 

-0.61 

-0.60 

As  may  be  seen  from  Table  35,  the  calculated  and  measured  values  com¬ 
pared  very  closely,  and  therefore,  confirm  the  accuracy  of  the  transfer  function. 

Let  us  now  consider  the  case  when  is  off  and  Q2  is  non- saturated. 
This  condition  is  represented  by  the  DC  equivalent  circuit  presented  in  Figure  65. 
The  analysis  of  this  circuit  is  quite  sUnple  since  most  of  the  nodal  vbUages  are 
known,  i.  e. ,  the  vol^ges  V^|  ,  Vg|  ,  known  and  are  ehown  on  tins 


lit 


diagram.  Also,  Ig  is  approximately  equal  to  and  since  it  is  of  the  order  of 
10"^  amps,  it  is  ne?,lected.  Therefore,  the  analysis  yields  only  equation  (184), 


♦  12V  •♦12V 


Figure  65.  DC  Equivalent  Circuit  of  Flip-Flop  with  off  and  Qj  Non-Satural»d 


(184) 


Since  the  value  of  Ig  ii  negligible,  Vg^  may  be  solved  for  directly  and 
is  -0.  605  volts.  It  should  be  noted  that  the  voltage  Vg2  for  Q2  non- saturated 
is  the  same  as  Vg^  for  Q2  saturated,  thereby  pointing  out  that  Q2  does  not 
operate  in  the  non- saturated  mode  but  is,  in  fact,  cut  off. 

3.  3.  2  Transfer  Functions  of  the  SteerinK  Network 

The  circuit  analyzed  in  this  subsection  is  that  of  a  steering  and  logic  net¬ 
work  designed  to  provide  the  input  pulse  for  the  Microminiature  Flip-Flop  circuit. 
Each  wafer  contains  two  independent  RC  differentiating  networks  with  level  block¬ 
ing,  diode-gated  outputs.  Two  i  iput  caj>acitors  are  provided  on  each  gate  to  allow 
a  wider  variety  of  input  pulses  to  be  used.  The  schematic  of  this  circuit  is  pre¬ 
sented  in  Figure  66. 


4.SV 


FIgur*  66.  Micromlniatur*  Steering  Network  Schematic 

As  may  be  seen  from  the  schematic,  the  circuit  is  of  symmetrical  nature, 
and  therefore,  an  analysis  of  one  side  of  the  circuit  is  sufficient.  The  analysis  of 
this  circuit  was  made  with  consideration  of  the  previous  stage  (e.  g. ,  in  the  system 
represented  in  Figure  6  the  steering  network  is  fed  by  a  NOR  circuit).  Therefore, 
the  appropriate  circuitry  of  the  NOR  circuit  is  included  into  the  analysis  of  the 
steering  network.  This  approach  requires  that  the  steering  network  be  analysed 
when  the  NOR  circuit  is  in  the  saturated  and  non- saturated  modes.  This  must  be 
done  since  different  portions  of  the  NOR  circuit  in  its  different  modes  have  a  signi¬ 
ficant  effect  on  the  steering  network. 

We  will  first  consider  *:he  case  when  the  NOR  circuit  (presented  in  Figure 
69)  is  operating  in  the  saturated  state.  The  DC  equivalent  circuit  for  this  condition 
is  shown  In  Figure  67. 

Referring  to  Figure  67(a),  R^  is  the  resistance  of  the  diode  CR^  and  R^' 
is  the  collector  resistance  of  the  NOR  circuit.  Figure  67(b)  is  the  simplified 
circuit  where 

R  R  ' 

*^p  “  Rj+'R^' 

An  analysis  of  this  circuit  yields  equations  (186)  and  (187).  Since  there  are 
only  two  unknowns,  it  is  much  easier  to  solve  these  equations 


120 


Figure  67.  DC  Equivalent  Circuit  of  Steering  Network  when  NOR  Circuit  Is  Saturated 


(186) 


(187) 


by  simple  substitution  rather  than  by  the  use  of  a  computer.  This  type  of  approach 
gives  the  following  relationships: 


V,(R.R.  +  R  R.  +  R  R. )  -  12  R  .R^ 
l'A5  p5  pA'  A5 

- - 

A  p 


V 


1 


s 


+ 


(188) 


(189) 


By  substituting  (189)  into  (188)  the  following  es^res Sion  for  as  a  fui^ticm  the 
circuit  components  is  derived: 


(190) 


^0  = 


4.  5  (Ry^Rj  +  RpRs  +  *^3*^A*^ 


A  5  p5  pA  3A  3p 


Once  the  expression  for  V^  is  obtained  the  solution  for  is  also  known  by  sub¬ 
stituting  the  expression  for  into  equation  (189).  The  numerical  solution  is 
then  obtained  by  substituting  the  nominal  values  of  the  components  into  the  expres¬ 
sions.  This  procedure  yields  the  following  results: 


Vq  »  6.  34  volts 
Vj  =  5.  99  .volts 

Let  us  now  consider  the  case  where  the  NOR  circuit  is  operating  in  the  non« 
saturated  mode.  For  this  condition  the  equivalent  circuit  presented  in  Figure  68 
applies. 


Figura  66.  DC  Equlvsiant  CtroiH  of  Sloorhg  Notweik  whon  NCR  CInuII  b  Non-Scrtunriod 


Referrii^  to  tiio  oquivaloat  circuit,  tho  ventage  Vj  is  clssnpod  ot  S.  I  v^ts 
assuming  a  0. 6-volt  hop  across  tho  diodo  CR^.  It  is  now  possiUo  to 
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deacribe  the  voltage  in  terme  of  the  circuit  componenta  by  equation  (191). 
Simplifying  the  expreaaion 


0 


(191) 


and  aolving  for 


V 


0 


a 


Vq  yielda  the  following  relationahip: 
4.  5  Rg  +  5.  1  Rj 


(192) 


Subatltuting  the  nominal  values  for  the  componenta  into  the  expreaaion  yielda  the 
voltage  Vq  when  the  NOR  circuit  is  operating  in  the  non- saturated  mode. 


Vq  =  5.0  V 

Thus,  we  have  completely  described  the  steering  network  in  terms  of  its  transfer 
functions  for  both  internal  and  external  effects. 


3.  3.  3  Transfer  Functions  of  the  Microminiature  NOR  Circuit 


The  circuit  to  be  analyzed  in  this  sub- section  is  that  of  a  NOR  circuit. 

This  circuit  is  simply  a  gate  with  a  logic  sense  of  plug  5  volts  »  False  and  zero 
volts  -  True.  The  circuit  is  composed  of  a  three-input  "OR"  gate  and  an  invertixm 
transistor  amplifier  and  a  clamping  diode  for  greater  reliability.  The  circuit 
accepts  a  positive  pulse  at  the  input  and  delivers  a  negative  pulse  to  the  succeeding 
stage.  The  NOR  circuit  schematic  is  presented  in  Figure  69.  The  performance 
criteria  of  the  circuit  is  shown  in  Table  36. 

This  circuit  was  analyzed  for  the  transistor  operating  in  both  the  saturated 
and  non- saturated  modes. 

We  will  now  consider  the  circuit  when  the  transistor  is  operating  in  the 
non- saturated  mode.  The  DC  equivalent  circuit  for  this  condition  is  preseitfed  in 
Figure  70. 


TABLE  36 

PBtFORMANCE  CRITRIA  OF  MICROMINIATURE  NOR  CKCUIT 
INPUT  REQUIREMENTS! 

Meodmum  "Fan  tn*^*  6  (oddttien  of  3  dledM  to  pint  5  or  12) 

Minimum  input  puiio  width  ■  40  m 
"faiao"  voitago  iovoi  ■  *4,5  to  +  6.0  volti 
"truo"  voitago  Iovoi  *  0  to  *0,5  volh 

OUTPUT  REQUIREMENTSt 

R|_  min.  (Max.  load  to  ground)  ■  0.5K  (1  mo) 

RS  min.  (Max.  load  to  ^CC)  *  ^  ^ 

Maximum  DC  "fan  out”  *  3  logic  inpuh  or  2  iovoi  Inputa 

Iovoi  goto  37AA02) 

*  Mox.  copoeitlvo  bod  *  100  pf 

Max.  AC  "Fan  out"  ■  4  puiao  inputa  (pulao  Iovoi  goto  37AA02) 

Typical  output  charactoriattca  for  oporation  with  5  volta  input 
puiao  at  25^  unioadod 
Turn  on  tbno  *  20  na 
Turn  off  tbno  *  20  na 
Turn  on  doloy  ■  10  na 
Turn  off  doi^  «  50  na 

In  order  for  the  traniistor  to  be  physically  in  non- saturation,  one  of  the  input 
diodes  has  to  be  grounded.  Referring  to  the  equivalent  circuit,  the  analysis 
yields  equations  (193)  through  (195)« 


V 


c 


5.  1 


-  4  -  V, 


B 


12  -  V 


1 


(193) 

(194) 

(19i) 


A 

However,  in  equation  (194),  Ig  ^  ~  10  thus,  since  the  effects  ^  uro 

negligible. the  equation  may  be  rewritten  as  equation  (196). 


4+  V 


B 


3 


(lli) 


Utilising  the  above  equations  the  matrix  equation  (197)  may  be  formed. 


1  1  1 

4 

^1 

S 

111  1 

qr 

12 

^  ^  ^ 

(197) 


The  solution  of  this  matrix  along  with  the  values  of  the  circuit  parameters  meas¬ 
ured  on  a  breadboard  in  the  laboratory  are  presented  in  Table  37. 


TABLE  37 

COMPARISON  OF  MEASURED  AND  CALCULATED  VALUES  OF  THE  NOR  CIRCUIT 
WITH  THE  TRANSISTOR  IN  THE  NON-SATURATED  MODE 


n 

n 

n 

mm 

■■ 

MecMured 

5,0 

0.8 

.... 

Calculated 

5.1 

0.73 

BSrI 

As  may  be  seeh  from  the  table,  the  calculated  and  measured  values  of  the  circuit 
compare  quite  closely  and  thereby  verify  the  accuracy  of  the  transfer  function. 


We  will  now  consider  the  case  where  the  transistor  is  operating  in  the 
saturated  mode.  The  DC  equivalent  circuit  for  the  condition  is  presented  in 
Figure  71.  Utilizing  this  circuit,  we  may  write  equations  (198)  through  (202). 


12  -  V, 


V,  -  V„ 


The  matrix  equation  ia  then  aoived  by  meana  of  a  computer  and  the  meaaured  and 
calculated  valuea  are  preaented  in  Table  38. 

TABLE  38 

COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF  THE  NOR  CIRCUIT 
WITH  THE  TRANSISTOR  IN  SATURATION 


An  examination  of  Table  38  reveala  the  cloae  comparison  between  the 
meaaured  and  calculated  valuea  of  the  circuit  parameters  thereby  confirming  the 
accxuracy  of  the  transfer  function  of  the  microminiature  NOR  circuit. 


SECTION  IV 


SYSTEMS  APPLICATIONS 


The  procedures  and  applications  described  in  this  section  are,  to  a  certain 
extent,  similar  to  those  discussed  in  Section  3.  The  major  differences  in  the 
technique  as  applied  to  systems  are:  (1)  the  added  complexity  and  (2)  the  prepara¬ 
tory  work  that  must  be  accomplished  in  order  for  the  results  to  be  both  meaning¬ 
ful  and  useful. 

The  study  has  yielded  three  possible  approaches  that  may  be  used,  all  of 
which  present  very  useful  results.  The  particular  approach  to  be  chosen  depends 
upon  a  number  of  factors  among  which  are:  (1)  type  of  information  required, 

(2)  complexity,  (3)  resources,  (4)  time  and  money  available,  etc.  Therefore, 
the  three  approaches  will  be  presented  along  with  a  discussion  of  the  afore¬ 
mentioned  factors.  The  first  method  that  will  be  discussed  is  the  most  elaborate 
and  consequently  the  most  expensive.  The  second  method  is  a  modification  of  the 
first  and  the  third  and  least  e3q>ensive  is  a  relatively  simple  but  quite  useful 
technique.  The  three  techniques  are  described  below. 

(1)  In  the  first  approach,  one  obtains  a  complete  transfer  function  in 
terms  of  all  the  components  in  the  system.  This  transfer  function 
may  be  derived  in  a  number  of  ways.  Two  methods  are  described 
below: 

(a)  The  transfer  fimctions  of  each  of  the  subsystems  may  be  individ¬ 
ually  obtained  and  then  by  adding,  subtracting  or  multiplyit^  Uie 
the  individual  transfer  functions,  depending  i;^n  the  function  of 
the  subsystems,  the  complete  transfer  function  may  be  obtained. 
However,  if  these  individxial  transfer  functions  are  in  matriac 
form,  the  rules  for  manipulating  matrices  must  be  followed,  i.  e. , 
in  order  to  multiply  matrices,  they  must  be  of  the  same  order, 
etc.  It  is,  therefore,  quite  probable  that  a  great  deal  of  prepara¬ 
tory  work  will  be  required  vnpon  the  individual  transfer  fanetimis 
prior  to  the  task  of  formulating  the  over-all  function.  This 
method  is  considered  to  be  the  easier  and  most  usdhal  method. 
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(b)  Another  method  that  may  be  used  to  determine  the  system  transfer 
function  is  to  eaqpress  it  analytically  in  terms  of  the  input  and 
output  voltages  of  the  system.  However,  with  this  type  of  ap¬ 
proach  a  great  deal  of  information  may  be  lost  since  many  of  the 
component  parameter  relationships  internal  to  the  input-output 
sections  of  the  system  may  not  appear  in  the  transfer  equation. 

It  is,  therefore,  recommended  that  the  first  technique  be  used. 

Once  the  system  transfer  function  is  obtained  and  is  formed  into  a 
transformation  matrix  of  the  order  N  ,  the  Monte  Carlo  process 
described  in  Section  II  of  this  report  may  be  used.  For  this  particular 
technique  one  would  run  the  process  for  every  component  and  combina¬ 
tion  of  components  in  the  matrix. 

(2)  The  next  approach  that  may  be  used  is  identical  to  that  described 
above  except  that  the  Monte  Carlo  process  is  applied  only  to  the 
critical  components  in  the  system.  An  example  of  this  technique  is 
described  below: 

Consider  the  simple  amplifier  shown  in  Figure  72.  It  is  obvious  from 
inspection  of  the  circuit  that  the  gain  is  given  by  Equation  (204); 


•  4 
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(204) 


Gain  ■  -2- 
*ln 


*m  ^ 


If  g  »  1  It  Is  easily  seen  that  equation  (206)  may  be  re- 
writfen  as  shown  below: 


Gain 


(205) 


It  may  be  stated,  therefore,  that  the  critical  components  of  the  ampli¬ 
fier  are  the  plate  and  cathode  resistances  R^^^  and  R|^  .  The  trans- 
conductance  may  be  neglected  since  it  is  common  to  both  numera¬ 
tor  and  denominator  of  the  gain  e:q>ression.  Therefore,  if  one 
determined  the  critical  components  of  a  system  it  would  be  possible  to 
save  a  great  dead  of  time  and  expense  using  this  approach  as  opposed 
to  that  described  in  (1). 


(3)  The  final  technique  to  be  described  is  one  in  which  the  Monte  Carlo 
process  is  applied  to  the  individual  subsystem  transfer  functions. 
With  this  type  of  approach  it  is,  therefore,  necessary  to  know  the  in¬ 
put  and  output  requirements  of  the  individual  subsystems.  Once 
the  requirements  are  known,  the  Monte  Carlo  process  may  be 
applied.  The  technique  consists  of  Monte  Carloing  each  of  the 
components  in  the  individual  subsystem,  and  noting  the  effects 
\q>on  the  output.  When  and  if  the  ou^ut  does  not  meet  the  ii^ut 
requirements  of  proceeding  subsystem,  the  values  of  the  com¬ 
ponent  or  components  are  noted  and  the  process  is  halted  in  the 
presence  of  a  failure  or  continued  on  the  next  subsystem  if  no 
fialure  occurs.  Another  useful  outcome  of  this  technique  is  that 
although  the  output  may  just  fall  within  the  input  req\iirements 
of  the  proceeding  stage,  it  may  cause  a  failure  at  another  point 
in  the  system  under  consideration.  As  an  illustration  ci  tUs 
technique,  let  us  consider  the  following  example: 


Suppose  we  have  two  amplifiers  in  series 


Flgwt73.  lUiNirariw  (EMMiftie 


The  gein  of  Amplifier  f  1  *  Amplifier  #2 
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The  Amplifier!  have  the  following  Requirement!: 


TABLE  39 

OPERATING  PARAMETERS  OF 
ILLUSTRATIVE  EXAMPLE 


Input 

Output 

Anpitner'l 
Anpilftsr  *2 

2^5%votH 
10+ 5%  volt* 

10+5%veih 

50  ^  5%  volts 

We  will  aesume  that  the  input  to  Amplifier  #1  ie  alwaye  at  2  volte.  The 
next  etep  ie  to  perform  a  Monte  Carlo  proceee  on  each  of  the  cotn^tm- 
ent!  in  the  amplifier  and  note  the  effects  upon  the  output.  From  the 
cumulative  distribution  functions  obtained  by  this  procedure,  it  is 
possible  to  determine  when  the  output  falls  outside  the  iiq^t  require- 
ments  of  the  proceeding  amplifier  and.  thus,  the  reliability  of  the 
amplifier  is  known.  This  process  is  then  carried  on  until  all  of 
the  subsystems  have  been  analysed. 


This  technique  is  considered  to  be  the  least  eiqMiasive  and  requires 
less  computer  time,  but  as  in  the  two  techniques  previously  men¬ 
tioned.  quite  time  consuming.  This  procedure  may  be  used  mi  all 
components  or  only  on  the  critical  components.  It  is,  therefore, 
recommended  that  this  technique  be  used  since  it  provides  all 
necessary  informsdion. 

The  techniques  outlined  above  will  provide  the  followiflf  iaformatioa: 

(1)  Based  on  the  underlying  frequency  distribution  ai  tihe  components, 
it  is  possible  to  dalerttttee  tte  reliability  of  a  systma  at  aero  Uwm 
a*  Mime  Ume  mr  ttrnes  in  future. 


(2)  It  provides  a  means  of  improving  performance,  e.  g. ,  let  ua  etmsider 
a  circuit  whose 


(«)  (b) 


Figure  74.  Illwtratlve  Example  of  Application  of  Monte  Carlo  Technique 

performance  is  plotted  as  the  cumulative  distribution  function  of 
Figure  74b.  The  normal  distribution  shown  in  Figure  74a  represents 
the  critical  conq>onent  R  of  the  circuit.  The  particular  tolerance 
limits  on  this  con^nent  (4  Tj)  will  yield  a  reliability  of  1  •  (X  4  X|) 
for  the  circuit.  However,  it  is  easily  seen  that  if  a  stricter  toler* 
anee  (Tg)  is  imposed  on  the  component  R  the  performance  of  Uie 
circuit  may  be  improved  as  represented  by  the  dashed  curve  in 
Figure  74b.  Therefore,  by  utilising  the  technique  in  this  manner, 
the  over -all  performance  of  a  circuit  or  system  may  be  inq^roved. 

(3) .  With  the  information  obtained  from  this  technique,  it  will  be 

possible  to  plan  better  spare  parts  allocation  and  maintainability 
methods. 

(4)  It  will  provide  a  better  method  for  determinii^;  trouble -sboethi^ 
techniques  since  the  critical  parts  will  be  highUghted  nd  tiM  system 
degradation  due  to  the  failure  of  tiisss  parts  will  be  known. 


4. 1  TIMING  SYSTEM 


It  is  the  purpose  of  this  subsection  to  present  a  specific  application  of  the 
Monte  Carlo  Technique,  previously  described  in  Section  IV,  on  a  system.  The 
system,  as  presented  in  Figure  5,  generates  a  constant  10  microsecond  pulse  at 
a  repetition  rate  of  100  kc  in  an  environment  where  a  large  shift  of  nominal  part 
parameters  may  be  experienced.  This  treatment  presents  the  specific  application 
of  the  technique,  however,  in  keeping  with  the  purpose  of  this  report,  that  of 
developing  a  new  technique,  the  routine  mechanics  are  not  included.  The  mathe¬ 
matical  simulation  procedure  is  as  follows: 

An  oscillator  is  synthetically  constructed  utilizing  the  transfer  func¬ 
tion  previously  developed.  The  critical  performance  criteria  of  the 
oscillator  is  the  resonant  frequency  and  the  satisfaction  of  this  cri¬ 
teria  will  automatically  ensure  adequate  input  to  the  high  level  am¬ 
plifier.  Once  it  has  been  ascertained  that  the  oscillator  is  operating 
within  specifications,  the  second  stage  or  high  level  amplifier  is 
synthetically  constructed  and  the  process  is  repeated.  This  proce¬ 
dure  is  followed  with  each  of  the  remaining  stages  in  the  system. 

If  a  failure  is  experienced,  the  process  is  immediately  halted  and 
the  cause  and  location  of  the  failure  noted.  The  whole  process  is 
repeated  1000  times  thereby  yielding  the  reliability  of  the  system. 
Referring  to  the  test  system  and  the  waveform  development  shown 
in  Figure  5,  the  following  requirements  must  be  satisfied  by  the 
respective  subsystems  in  order  to  obtain  successful  operation. 
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4.  Z  0B8CIUPT10N  OF  THE  M1CROII1NIA.TUIIE  ORCUXTS  IN  THE  ERRCA 
SENSING  AND  READOUT  CIRCUIT  '  ... 
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In  Section  3. 3  tho  throo  mierominiaturo  circuits  nrs  doscribsd  which 
compose  the  systems  in  Section  IV.  and  detailed  schematic  diagrams  are  drawn. 
In  this  subsection  these  same  three  circuits  are  again  described,  but  hers  in 
terms  of  their  logical  functions  or  o£  their  general  function  in  a  configuration, 
as  opposed  to  their  electrical  configuration.  [26,  27l 

4.  2. 1  Logical  Description  <rf  a  Circuit 

All  electrical  circuits  can  be  divided  into  two  types,  digital  and  analog. 

If  a  circuit  is  viewed  as  a  black  boa,  i.  Sm  only  tbs  inputs  and  ou^nits  are  of  ' 
interest,  then  analog  circuits  are  those  in  which  the  voltages  at  the  terminals 
may  vary  continuously  over  a  range;  whereas  digital  circuits  are  those  in  which 
the  terminal  voltages,  both  input  and  ou^ut,  are  restricted  to  two  (or  any 
finite  number  but  usually  only  two)  sufficiently  separated  narrow  sub-ranges, 
say  Rj  and  R2*  of  the  total  possible  range  R.  If  the  voltage  falls  outside  of  R^ 
and  R2  the  circuit  is  said  to  have  failed.  If  the  circuit  dose  not  fail,  then  the 
voltage  at  each  terminal  must  be  in  one  and  only  one  of  two  mutually  esclusive 
sub-ranges.  Different  letters  (  a,  b,  — )  of  ti>e  alphabet  may  be  assigned  to 
each  terminal  to  indicate  the  presence  of  a  voltage  in  one  of  the  sub-ranges,  and 
the  negation,  in  the  senes  of  syinbolic  logic,  or  Boolean  algebra  of  this  letter 
(A,  B,  — )  means  that  the  voltage  is  in  the  other  sub-range.  An  mnuoaple  of  this 
is  given  in  Figure  7  5. 
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In  Figure  75  the  Symbolic  Logic  Block  Diagram,  the  following  voltagee 

define  ABC  and  A  B*  ^ 

It^ut  ■  A,  when  0  <  <  1 

Iiq;>ut  ■  A,  when  5  <  <  6 

li^ut  ■  B,  when  0  <  <  1 

Input  ■  B|  when  5  <  <  6 

Ii^ut  ■  C(  when  0  <  V  <  1 

c  •• 

lapnt  •  C,  when  5  <  <  6 

,  A  non- sequential  circuit  is  one  in  wl^ch  the  outputs  depend  only  on  the 

inputs  and  nothing  else,  e.gi,  internal  state,  past  history,  etc.  In  a  non¬ 
sequential  circuit  the  same  inputs  always  yield  the  same  outputs.  A  sequential 
circuit  is  one  in  which  the  outputs  depend  not  only  on  the  inputs  but  also  on  the 
internal  state.  Different  internal  states  handle  the  same  iiq>uts  differently  i.  e., 
may  yield  different  outputs  for  the  same  set  of  inputs.  In  this  sense  the  in¬ 
ternal  state  may  profitably  be  regarded  as  itself  another  input.  Sequential 
circuits  are  sometimes  described  as  those  having  "memory",  since  they 
"remember"  past  inputs  by  correlating  different  internal  states  with  various 
sequences  of  past  inputs.  A  complete  description  of  a  circtiit  can  be  given  by  a 
truth  table,  if  the  inputs  and  outputs  are  represented  by  their  respective  letter 
symbols  and  arranged  in  the  following  manner.  The  letters  are  listed  horiacm- 
tially  as  column  headings  and  each  row  gives  each  possible  state  of  the  circuit  . 
Affirmation  and  negation  can  be  described  as  in  symbolic  logic  (Boolean  A^ebra) 
or  numbers  can  be  used  (  'one*  for  affirmation,  'aero'  for  negation).  The 
resulting  array  is  called  a  truth  table  (for  example  Table  40)  and  is  a  ccnnplete 
description  of  all  possible  states  of  a  non- sequential  digital  circuit.  A  Boolean 
equation,  or  a  logical  truth  function,  is  completely  equivalent  to  this  array 
and  represents  the  fact  that  the  ouq^s  of  the  circuit  are  a  certain  function  of 
the  ii^ts.  To  say  that  a  digital  circuit  performs  a  certain  lexical  operation 
or  that  it  has  a  certain  logical  descriptUm  is  the  same  as  saying  tiiat  Its  inputs 
and  ouq^s  satisfy  a  particular  Boolean  equation.  itoMe  the  name  of  sqn^ 
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ia  tha  same  aa  saying  that  its  inputs  and  outputs  satisfy  a  particular  Booleui 

0 

equation.  Hence  the  name  of  some  digital  circuits  (e.  g.  NOR  circuit,  NAND 
circuit)  often  comes  from  the  kind  of  equation  which  describes  its  operation. 

(e.  g.  (A  B)  for  not  (A  or  B)  ) 

4.  2.  2.  NOR  Circuit 

The  NOR  circuit  used  here  has  3  ii^uts  A,  B,  C  and  one  output  D.  The 
truth  table  below  is  arranged  according  to  the  explanation  given  above.  For  all 
possible  input  configurations  the  output  D  is  given  as  a  function  of  this  con¬ 
figuration  (  of  ordered  triplets  of  'aero*  and  'one'). 


TABLE  40 

REPRESENTATIVE  TRUTH  TABLE 


A 

B 

c 

D 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

1 

0 

1 

0 

0 

1 

1 

0 

1 

1 

1 

0 

The  equation  equivalent  to  this  truth  table  is  given  by  equation  (206L 

(A  +  B  +  C)  ■  A  B  C  »  D  (204) 

4.  2. 3  Flip«Flop  Trigger  Circxtits 

The  conibination  of  the  flip  flop  and  trigger  circuit  is  dctdgBCd  to  iHt 
as  a  single  unit  and  forms  a  sequential  circuit. 
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The  flip  flop  performs  an  operation  which  is  often  referred  to  as 
"remembering".  It  has  two  possible  complementary  outputs,  called  'true*  and 
'false'.  It  can  be  switched  from  one  state  to  another  by  short -duration  pulses 
and  "remembers"  indefinitely  the  last  state  into  which  it  was  thrown.  If  the 
flip  flop  is  in  the  condition  where  the  'true'  output  is  in  the  'aero'  state  and  a 
pulse  input  is  delivered  to  the  'true'  side,  the  flip  flop  will  change  state  so  that 
the  'true'  side  goes  from  'one'  to  'zero'.  Consequently  the  'false'  side  goes 
from  'zero'  to  'one'.  A  signal  last  delivered  to  the  'false'  input  will  settle  the 
flip  flop  so  that  the  'one'  side  goes  from  'one*  to  'zero',  and  the  'zero'  side 
generates  a  'true'  signal  (i.  e. ,  goes  from  'zero'  to  'one').  This  state  is  main¬ 
tained  indefinitely  and  is  referred  to  as  the  'one'  state  of  the  flip  flop.  A  second 
short -duration  pulse  delivered  to  the  'true'  input  of  the  flip  flop  does  not  change 
the  state,  but  a  pulse  delivered  to  the  'false'  input  reverses  the  state  of  the  flip 
flop  so  that  its  'true'  output  terminal  generates  a  'zero'  signal  amd  its  'false' 
side  a  '  one'  signal.  The  flip  flop  is  now  in  the  '  zero'  state  and  remains  there 
until  an  input  is  next  delivered  to  the  'false'  input  terminal.  See  figure  76. 


Figure  76.  Logic  Dtogram  of  Flip-Flop 

The  trigger  circuit  (see  Figure  77)  is  composed  of  two  identical  circuits. 


Figure  77,  Logic  Dfc^pom  of  Tr^gor  Clreuit 
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Its  two  outputs  ars  connaeted  to  the  two  inputs  of  ths  fl^-fU^  circuit.  THs 
system  clock  is  coansetsd  to  ths  ii^uts  of  ths  tr^sr  circuit.  Clock  pulsss  pass 
from  ths  trigger  input  to  ths  trigger  output  terminal  if  ths  level  ii^put  is  in  stats 
'one'.  If  ths  level  input  is  in  state  'aero',  ths  clock  pulses  do  not  pass  through  to 
ths  output. 

4.2.4  Description  of  ths  System 

The  function  of  ths  error  sensing  and  readout  system  it  to  monitor  two 
inputs,  each  input  being  an  ou^ut  from  one  of  two  identical  systems,  utd  to  com* 
pare  these  ii^uts  with  a  reference.  An  indication  signal  is  provided  when 

1.  Either  one  of  the  two  systems  is  in  error  (but  only  one)  or 

2.  Both  systems  are  in  error. 

The  error  sensing  subsystem  is  a  gating  structure  composed  entirely  of 
NOR  circuits  which  has  one  indicator  for  single  error,  and  one  indicator  for 
double  error,  so  that  correct  operation  of  the  monitored  systems  corresponds  to 
an  absence  of  such  indications.  Correct  operation  is  defined  as  agreement  with  the 
reference  signal  at  all  times. 

The  readout  subsystem  is  composed  of  two  flip  flop  circuits  that 

"remember"  that  an  error  has  occurred  in  the  time  interval  following  the  last 

reset  of  the  system.  These  flip  flops  provide  an  indication  signal  at  the  time  of 

malfunction  and  continue  to  supply  this  indication  until  the  system  is  reset.  Re- 

setting  deers  the  system  of  all  error  indications.  Since  these  are  digital  systems* 

the  signals  at  any  time  will  be  two  valued  (as  has  been  described  above).  Ther^oro 

if  the  two  systems  being  monitored  have  different  ou^uts  at  any  time*  one  must 

agree  with  the  reference  and  the  other  must  disagree.  If  the  'true'  outputs  of 

_  • 

these  systems  are  represented  by  A  and  B  then  the  requirement  for  single  failure 
(S)  can  be  described  by  the  Boolean  equation 

S  ■  A  B  4  a'B  {7m) 

la  the  case  of  double  failure  (D)*  the  two  systems  muM  agree  u^  Mph 
other  and  disagree  with  the  reference  (R).  TIm  equattea  for  this  te 


(208) 
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D-ABR-I-ABR 

It  follows  that  correct  operation  (C)  means  that  the  reference  and  both 
systems  agree.  This  can  be  represented  by 

C>ABR  +  ABR  (209) 

The  equations  for  S  and  D  can  be  rewritten  in  the  following  manner 

S  *  A  B  +  A  B 
S  =  (A  B  +  A  B) 


S*(AB)  (AB)  »  [jA  +  B)  +  (A  +  B) J  (210) 

D  =  ABR  +  ABR+(AB  +  AB)(RA  +  RA) 

D  »  (A  B  +  A  B)  (R  +  A)  (R  +  A) 


D=  [(A  B  +  A  B)  +  (R  + A)  +  (R  +  A)] 

But  (A  B  +  A  B)  =  S 

(211) 

Therefore. 

D  »  [s  +  (A  +  R)  +  (A  +  R)J 

(212) 

E]q>ressed  in  this  form  these  equations  describe  the  error  sensing  NOR 
circuit  configuration  of  Figure  6. 

In  the  error  readout  section  (also  shown  in  Figure  6)  the  output  of  T^  will 
trip  flip  flop  F2  on  receiving  the  first  clock  pulse  during  the  time  Bone  in  which 
S  is  present.  In  a  similar  manner,  T^  will  trip  F^  on  receiving  the  first  clock 
pulse  occurring  while  O  is  present.  The  flip  flop  drives  the  appropriate  error 
indicator  when  S  or  O  have  occurred. 

The  Monte  Carlo  techniques  described  in  Section  4. 1  are  directly  ap¬ 
plicable  at  this  point. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  report  has  developed  mathematical  simulation  procedures  to 
accomplish  the  following: 

1.  Provide  the  designer  with  a  specified  degree  of  confidence  that  his 
circuit  will  have  a  particular  performance  criteria. 

2.  Provides  a  means  of  readily  evaluating  the  sensitivity  ol  the  circuit 
thereby  indicating  the  major  sources  of  variability. 

3.  Tolerances  for  part  parameters  or  circuit  elements  can  be 
realistically  specified  for  the  selection  of  parts. 

4.  The  system  is  optimized  for  a  given  cost. 

5.  The  development  of  transfer  functions  require  that  the  definition  of 
failure  be  precisely  specified  mathematically  and  thus  the  usual 
decisions  as  to  the  acceptability  of  a  circuit  performance  criteria 
is  eliminated  . 

6.  The  cost  of  the  parts  purchased  can  be  reduced  as  the  worst  case 
design  technique  is  no  longer  required  and  the  parts  do  not  have  to  be 
as  precise. 

7.  The  precision  of  this  technique  is  limited  only  by  the  number  <d 
iterations  performed;  however,  this  is  non-restrictive  because  of  the 
availability  of  modern  day  computers. 

8.  The  technique  makes  possible  the  planning  of  inqproved  allocation  of 
spare  parts. 

9.  It  is  possible  to  improve  maintainability  techniques;  i.  e. ,  the  infor¬ 
mation  will  be  available  to  allow  the  optimum  determination  of  repair 
crews  required  to  maintain  the  system. 

10.  It  will  provide  a  guide  for  the  determination  of  trouble  elMottag 
techniques  since  the  critical  parts  will  be  highlighted  and  ^  sy^mn 
degradation  due  to  the  failure  of  these  parts  will  be  kaotm. 


11.  Bated  i^on  the  underlying  frequency  diatribution  of  the  componente, 
it  it  pottible  to  determine  the  reliability  of  a  ay  item  at  xero  time 
and  at  tome  time  or  times  in  the  future. 

Thus,  it  it  concluded  that  in  mott  catet  Mathematical  Simulation  it  the 
only  featible  aolution  to  the  reliability  prediction  problem. 

Recommendation!  for  future  work  in  thete  areat  include: 

1.  A  hardware  program  be  continued  to  further  validate  the  predictive 
techniques  developed  in  this  report. 

2.  Additional  work  be  continued  in  developing  the  optimum  diatribution 
of  parts  and  components  for  optimising  system  reliability  against 
given  constraints. 

3.  Research  be  conducted  into  the  further  application  of  the  technique 
in  the  area  of  automatic  circuit  design. 


GLOSSARY 


The  set  of  definitions  presented  in  the  RADC  Reliability  Notebook,  dated 

30  October  1959  (see  reference  94  in  the  Bibliography)  is  appropriate  for  use  in 

understanding  this  report.  If  an  alternative  definition  is  needed  for  clarity,  the 

list  below  may  prove  useful. 

FREQUENCY  •  The  ratio  of  the  number  of  events  which  meet  a  performance 
requirement  to  the  total  number  of  events. 

FIRST  FAILURE  TIME  mSTRlBUTION  -  The  distribution  of  all  times  to  the 
first  failure. 

DENSITY  FUNCTION  -  The  first  derivative  of  the  distribution  function. 

PART  •>  The  smallest  basic  element  of  a  complete  system. 

MODULE  -  (Normally  used  interchangeably  with  component. )  A  module  is  an 
article  which  is  normally  a  combination  of  parts,  subassemblies,  or 
assemblies,  and  is  a  self-contained  element  of  a  complete  operating 
equipment,  and  performs  a  function  necessary  to  the  operation  of  that 
equipment. 

CUMULATIVE  DISTRIBUTION  FUNCTION  -  (See  Appendix  A) 

SYSTEM  PERFORMANCE  CRITERIA  -  A  system  performance  criterion  is  a  rule 
that  states  within  what  limits  a  system  performance  measure  must  fall. 

A  system  performance  measure  may  be  of  two  types,  a  deterministic 
performance  measure  or  a  probabilistic  performance  measure.  A  ^er- 
ministic  performance  measure  is  a  measure  that  can  be  characterised  by 
a  single  nxunber. 

Some  examples  of  this  type  of  measure  are: 


a. 

gain 

d. 

rise  time 

b. 

bandwidth 

e. 

fall  time 

c. 

signal  to  noise  ratio 

f. 

telay  time 

g.  noite  figure  1.  output  impedance 

h.  input  impedance 

A  probabilistic  performance  measure  is  a  measure  that  can  be 
characterized  by  a  cumiilative  distribution  function  which  may  take  any  of  the 
following  forms: 

1.  an  equation 

2.  a  graphical  curve 

3.  an  approximating  polynomial 

CRITERION*-  A  standard  of  judging,  a  rule  or  test  by  which  anything  is  tried 
in  forming  a  correct  judgment  respecting  it. 

SPECIFICATION*  -  Minute  description  of  particulars  or  the  particular  details 
themselves,  e.  g./in  the  terms  of  a  contract  the  details  of  construction 
etc.  According  to  the  glossary  in  the  RADC  Reliability  Notebook,  a  detailed 
description  of  the  characteristics  of  a  product  and  of  the  criteria  which 
must  be  used  to  determine  whether  the  product  is  in  conformity  with  the 
description. 

NOMINAL  VALUE  -  The  stated  value  (e.  g. ,  by  the  manufacturer  or  designer) 
of  some  characteristic  or  measure  of  performance  of  a  piece  of  equip¬ 
ment.  Note:  Due  to  internal  and  external  stresses,  it  is  seldom  the  case 
that  the  nominal  value  coincides  with  the  actual  measured  value  of  a 
particular  characteristic. 

*  Definition  from  Webster's  New  Collegiate  Dictionary. 


V 


BIBLIOGRAPHY 


(1)  Matrix  Analyla  of  Electric  Network*  by  P.  L«  CorbciUar,  Harvard 
University  Prase.  1950. 

(2)  A  Treatise  on  Electricity  and  Magnetism  by  J.  C.  Maxwell,  Oxford,  ed 
3.  1892. 

(3)  Mathematical  Simulation  for  Reliability  Prediction  by  I.  Bosinoff,  et  al. 
Final  Report,  30  October  1962  (153  p)  illus.  Report  No.  F491>1,  Contract 
No.  AF30(602)«2376.  Unclassified  Report. 

(4)  Modern  Mathematics  for  the  Enflineer,  E.F.  Beckenbach,  ed. ,  McGraw 
Hill,  1956.  See  esp.  Chapt.  1 2  for  Monte  Carlo  methods. 

(5)  Modern  Probability  Theory  and  its  Applications,  E.  Parsen,  Wiley,  1960. 
see  esp.  pp  220>1  for  table  of  probability  laws  from  which  random  num« 
bers  may  be  generated. 

(6)  Biometrika  Tables  for  Statisticians  by  E.  S.  Pearson  and  H.  O.  Hartley, 

Vol.  1,  Cambridge  University  Press,  1956. 

(7)  Pulse  and  Digital  Circuits  by  J.  MiUman  and  H.  Taub,  McGraw  Hill,  1956. 

(8)  Reliable  Preferred  Solid-State  Functional  Divisions,  Technical  Report, 
RADC-TR«59-243,  15  December  1959,  Job  2131,  KLX-10177.  Contract 
No.  AF30(602)-1906,  prepared  for  Rome  Air  Development  Center,  Air 
Research  and  Development  Command,  U.  S.  A.  F.  This  contains  a  list  of 
conventional  circuits. 

(9)  IBM  Machine  Analysis  of  Network  of  Its  Applications  by  F.  H.  Branin, 

Jr. ,  March  30,  1962,  TR  00.  855, 

(10)  A  C<Maaputer  Application  to  Reliable  Circuits  Design  by  L.  Hellermaa, 

1  TR  6o.  786,  IMvelopment  Laboratory,  Data  Systems  IHvieicm,  IBM, 

New  York. 

(11)  The  Nature  of  the  Monte  Carlo  Method  by  J.  Neilson,  Technical  r^ort 

137,  November  5,  1951,  Cr^  Laboratory,  Harvard  Ualveraity; 

Office  ol  Naval  Research,  Contract  N5  ORI-66  Task  Order  1  NR  *  678-01 1. 

(12)  Vacuum  Tube  Oscillators  by  Edson,  W.  A. ,  John  Wiley  and  Sons  Iim. ,  1953. 

(13)  NPN  Triple  Diffused  8ilic<»  Me iw  Transistors  ad/89,  PaeiRc8«nl- 
conductors,  Inc. ,  Culver  ^ty,  Callionda. 


m 


(14)  Fairchild  NPN  Diffused  Silicon  Trarsistors  SL-5/3.  May  1961,  Fairchild 
Semiconductor  Corporation,  Mountain  View,  California. 

(15)  Control  Syetem  Syntheeie  by  J.  G.  Truxal,  McGraw  Hill,  Page  35. 

(16)  Reference  Data  for  Radio  Engineers,  I  T  and  T  Corporation,  4th  Edition, 
c  1949,  Page  1090. 

(17)  The  Kolmogorov~Smirnov  Test  for  Goodneea  of  Fit,  F.  J.  Massey,  Jour¬ 
nal  of  the  American  Statistical  Association  No.  46  Pages  68-78,  c  1951. 

(18)  Nximerical  Tabulation  of  the  Distribution  of  Kolmogorov' s  Statistic  for 
Finite  Sample  Siae,  2.  W.  Bimbaum,  Journal  of  the  American  Statistical 
Association,  No.  47,  Pages  425-451,  c  1952. 

(19)  Advance  in  Reliability  Prediction  by  I.  Bosinoff,  R.  Jacobs  and  J.  Herman, 
1962  National  Space  and  Telemetry  Conference,  Miami,  Fla. ,  October 
2-4,  1962. 

(20)  Physics  of  Failure  in  Electronics,  Rome  Air  Development  Center  and 
Armour  Research  Foundation,  September  26-27,  1962. 

(21)  Wave  Generation  and  Shaping,  by  L.  Strauss,  McGraw  Hill  -  1960- 
Page  258. 

(22)  Linear  Signal-Flow  Graphs  and  i^iplication,  by  C.  Chow,  John  Wiley  and 
Sons,  Inc. ,  New  York  1962. 

(23)  Signal  Flow  Graphs,  by  R.  F.  Haskins,  M.  Sc.,  Electronic  and  Radio 
Engineers,  August  1959. 

(24)  Stabilization  of  Pulse  Direction  in  Monostable  Multivibrator  by  A.  Luther, 
Jr.  R.  C.  A.  Review,  1955. 

(25)  An  Analysis  of  the  Modes  of  Operation  of  a  Simple  Transistor  Oscillator 
by  J.F.  Gibbons,  Tech  report  1713-1,  Proceedings  of  the  1.  R.E. ,  Vol. 

49,  No.  9,  pp.  1383-1390,  Sept.  1961. 

(26)  Logical  Design  of  Digital  Computers,  M.  lister,  John  Wiley  and  Sons, 
1958. 

(27)  Design  of  Transistorised  Circuits  for  Digital  Computers,  A.L  Pressman, 
John  F.  Rider,  PubUeher,  Inc.  1959. 


APPENDIX  I 


DEFINITION  OF  CUMULATIVE  DISTRIBUTION  FUNCTION 

The  function  F(x)  ie  called  the  cumulative  distribution  function  (cdf)  of 
the  random  variable  X  if  and  only  if  F(x)  has  the  following  properties: 

(1)  F(x)  lies  between  0  and  1  for  all  values  of  x. 

(2)  F(x)  decreases  to  zero  as  x  decreases  to  -<o. 

(3)  F(x)  increases  to  1  as  the  upper  bound  b  of  the  interval  in  which 

X  is  supposed  to  lie  increases  to  +«> .  (The  probability  of  a  value  less 
than  +«  is  1. ) 

(4)  F(x)  is  an  increasing  function  in  the  sense  that ^f  a  b,  where  a  is 
the  lower  bound  of  the  interval,  then  F(a)  ■  F(b).  (The  probabi¬ 
lity  that  a  random  value  X  is  smaller  than  a  is  less  than  the 
probability  that  a  random  value  X  is  smaller  than  the  larger 
number  b. ) 

(5)  F(x)  is  always  continuous  on  the  left;  i.  e.  ,  F(a)  increases  to  F(b) 
as  a  increases  to  b. 

Conversely,  any  real  function  F(x)  with  the  above  properties  defines  a 
cumxilative  distribution  function  for  a  random  variable  X,  that  is 

F(x)  m  Pr(X  <  x) 

An  intuitive  grasp  of  a  cumulative  distribution  function  may  be  obtained 
from  the  following  illustration.  Siq>pose  a  thousand  resistors  have  been  made  and 
are  tested  to  find  the  exact  resistance  of  each  one.  If  the  thousand  values  resulting 
from  the  test  are  ordered  according  to  magnitude,  starting  with  the  smallest  and 
ending  with  the  largest,  a  graph  is  obtained  at.  shown  below  in  Figure  1-1. 

-  Here  the  ordinate  represents  the  number  of  resistors  having  a  value  less 
than  or  eqxial  to  that  shown  on  the  horixontal  scale.  Depending  on  the  scale,  ^e 
ordinate  is  sometimes  called  the  probability  (0  to  1)  or  relative  frequency  (0  to  1) 
or  percentage  (0  to  100),  but  in  any  case,  it  represents  the  value  of  the  etcnnlative 
distribution  function  F(34  . 
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The  derivative  of  the  cvuniilative  dietribution  function  F(x)  is  called  the 
denaity  function  (d.  f. )  'f(x).  Since  the  c.  d.  f.  is  ecaled  to  go  from  0  to  1,  it  must 
be  the  case  that  the  area  under  the  d.  f.lis  exactly  equal  to  1. 

The  items  discussed  above  can  be  extended  to  two  or  more  dimensions. 
When  a  process  simultaneously  produces  two  (or  more)  numbers  at  random  at 
each  trial,  the  collection  of  pairs  of  random  numbers  (X^,  X2)  can  be  considered. 
If  the  probability  is  known,  with  which  a  pair  of  random  numbers  (X^,X2)  lie  in 
the  quadrant  given  by  the  two  inequalities  X^  <  x^  and  X2  <  X2,  the  joint 
cumulative  distribution  function  can  be  defined  by 

F(xj,  X2)  ■  Pr  pCj  <  Xj  and  X2  <  X2} 
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APPENDIX  U 

DISCUSSION  OF  RELIABIUTY  MEASURES  AS  BASED 
ON  THE  CUMULATIVE  DISTRIBUTION  FUNCTION 

The  definition  of  reliability  that  has  received  wide  acceptance  in  the 
literature  la: 

The  probability  of  performing  a  specified  function  without  failure 

under  given  conditions  for  a  specified  period  of  time. 

This  can  be  expressed  mathematically  as: 

Reliability  ■  R(t)  «  e‘^*  (II- 1) 

where  X  «  failure  rate  i.  e. ,  reciprocal  of  mean  time 

between  failures 

t  3  duration  of  mission  or  task  , 

This  definition  of  Reliability  defines  R  as  a  probability  which  is  a  function  of 
time  (t)  and  of  internal  and  external  conditions,  since  \  depends  on  the  environ¬ 
ment  as  well  as  on  internal  stresses. 

In  this  equation,  the  expression  e'^^  gives  one  type  of  probability  function 
which  can  be  used  as  a  measure  of  the  reliability. 

There  may  be  other  types  of  probability  functions  which  may  be  used  as 
reliability  measures  in  addition  to  the  type  of  probability  function  described  above. 

The  following  analysis  shows  how  the  formulation  of  reliability  in  (II- 1)  is 
related  to  the  "Monte  Carlo"  procedures  described  in  Section  II.  This  analysis 
will  deal  first  with  a  system  having  only  one  type  of  unit.  To  express  the  above 
reliability  measure  R(t)  in  terms  of  a  cumulative  distribution  function  FqU) 
measuring  the  probability  of  failure  let: 

F^(t)  -  1  -  R^(tj  «  1  -  e"^"*) 

so  that  F.(t)  +  R_(t)  «  1 

n  n 

where  n  •  1,  2  ...  « 


(n-2) 
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As  t  approaches  <»  ,  approaches  the  steady  state  distri'  Mon  func¬ 

tion. 

F  is  the  particular  cumulative  distribution  function  which  describes  the 

n  ^  th  th 

distribution  of  the  various  intervals  of  time  between  the  n  “  and  (i.  <  1)  “ 

failure. 

Let  the  cumulative  distribution  function  of  a  circuit  perfurmaio  >  ■•i-  .>sure 
at  t  =  0  appear  as  shown  in  Figure  II- 1. 


Figure  IM 


Note  that  in  Figure  II- 1  each  of  the  cumulative  distribution  functions  repre¬ 
sents  the  same  set  of  systems.  A  circuit  is  considered  to  have  failed  whenever 
the  performance  measure  falls  outside  the  preassigned  limit  of  satisfactory  opera¬ 
tion.  Therefore,  any  circuit  or  system  whose  performance  measure  falls  to  the 
right  of  the  critical  performance  limit  has  failed  (as  shown  by  vertical  line  in  . 
Figure  II- 1).  From  Figure  II- 1  it  is  possible  to  construct  Figure  U-2  by  noting 
the  times  at  which  failure  occurs.  These  are  found  at  the  intersections 
critical  performance  limit  with  the  time  dependent  cumulative  distribution  func¬ 
tions.  The  time-to-first-failurc  distribution  Fj(t)  is  defined  as  1  -  R^(t). 

This  is  shown  in  Figure  11-3.  It  is  important  to  note  here  that  in  many  instances 
R|(t)  may  not  be  in  exponential  form,  since  it  is  a  direct  function  of  the  particular 
family  of  cumulative  distribution  functions  curves  shown  in  Figure  U-l. 


Figura  H-2 


FIgura  11-3 


The  fact  that  one  may  either  repair  or  replace  units  which  have  failed  will 
influence  the  development  of  the  mathematical  models  used  to  describe  the  relia¬ 
bility  measures  of  the  system.  Furthermore,  in  the  case  of  repair  there  will  be 
one  type  of  model,  while  in  the  case  of  replacement  another  type  of  model  will  be 
necessary. 

In  the  case  of  repair,  the  simulated  change  in  performance  criteria  is 
represented  as  shown  in  Figure  11-4.  To  determine  the  distribution  of  times  to 
first  failure  (i.  e. ,  F^)  and  the  distribution  of  times  between  the  consecutive  fail¬ 
ures  n  and  n+ 1  (i.  e. ,  F^)  a  critical  performance  limit  is  defined  (horisontal 
line  in  Figure  11-4.  The  time  to  first  failure  is  that  time  at  which  the  curve 
describing  the  performance  measures  in  time  first  crosses  the  line  representing 
the  critical  performance  limit. 
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If  these  times  to  first  failure  are  ordered  according  to  magnitude,  they 
yield  the  cumulative  distribution  fiinction  shown  in  Figure  II- 5. 


Figure  11-5 


Figure  11-4  shows  how  several  parte  come  to  the  time-of-first  failure. 
Figure  II-6  takes  one  of  these  parts  and  shows  the  performance  of  the  part  passing 
the  critical  performance  limit  at  t^  (time-to-first-failure),  being  repaired, 
rising  in  performance  again  until  it  passes  the  critical  performance  limit  at  t2 
(time-to-second-failure),  and  so  on.  Figure  11-7  does  this  for  another  part. 

If  a  whole  aeries  of  such  graphs  are  imagined  for  each  part,  then  a  cumulative 
distribution  function  may  be  found  for  each  At^  C(At^  *  ^i  ”  ^^i  "  1.  e.  ,  each 

incremental  time  between  failures.  This  gives  F.  (the  cumulative  distribution 
function  for  time  from  the  (i  -  1)  to  the  i  -  th  failure  for  the  type  of  part  being 
tested  on  the  basis  of  the  performance  of  a  large  number  of  actual  tested  samples. 
Notice  that  the  time  required  to  repair  a  part  is  taken  as  aero  on  the  graphs. 
Although  this  is  rare  in  practice,  this  time  for  repair  does  not  influence  the 
distributions  in  question. 

Interfailure  Time  Distribution  with  Monitoring  Replacement  Policy 

In  the  case  where  the  system  is  monitored  and  failed  parts  are  replaced, 
a  second  mathematical  model  is  necessary.  In  Figure  II-8  the  solid  curves  repre¬ 
sent  the  cumulative  distribution  function  of  the  original  ptqralatio^  oi  parts,  wldch 
were  subject  only  to  degradatim  in  time  and  no  monitoring.  If  the  defective  fNarte 
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ar«  raplaead  at  intarvala  whan  tha  cumulativa  diatribution  function  wlU  ba  daacribad 

by  tha  dottad  curvaa,  i.  a. .  (1  -  X)  of  tha  unlta  hava  failad  at  tima  t  ■  t  . 

n 

(1  -  X)  of  tha  unita  hava  failad  at  t*1000  hours;  thaaa  occur  bacauaa  a 
critical  performance  limit  is  exceeded.  The  replacement  of  (1  -  X)  of  the  units 
taken  from  the  population  represented  by  the  cumxilative  distributive  fvmction  will 
result  in  the  cumulative  diatribution  function  represented  by  the  dashed  curve  (1) 
at  t  e  0.  This  replacement  would,  at  t  «  2000  hours,  cause  this  curve  to  shift  and 
be  represented  by  the  dashed  curve  (2).  This  curve  would  be  composed  of  (1  -  X) 
of  the  samples  from  the  distribution  at  t  »  1000  and  X^  of  the  samples  from  the 
distribution  at  2000  hours.  The  net  result  would  be  that  the  reliability  of  the 
samples  of  the  curve  (2)  would  be  increased.  The  increase  in  reliability  for  any 
point  on  the  curve  would  be  a  maximiun  of  ( 1  -  X).  This  procedure  may  be 
repeated  indefinitely  as  time  increases. 


Figure  11-9 

Figure  Il>9  represents  a  plot  of  the  Reliability  as  a  function  of  time  with 
and  without  monitoring.  Figure  II- 10  is  the  corresponding  plot  o£  the  failure  rates. 

Figures  11-8,  11-9,  11-10  are  related  to  each  other  in  the  same  way  as 
Figures  U-1,  U-2,  U-3. 


System  with  Two  Types  of  Units 

To  obtain  the  over-all  system  reliability  for  a  system  with  two  typeis  of 
unita  it  is  sufficient  to  multiply  the  two  reliabilities. 


0-6 


Over-all  system  failure  occurs  If  any  one  of  the  two  units  fails.  Upon 
failure,  it  is  assumed  that  the  failed  unit  will  be  restored  to  an  operable  condition. 
For  simplicity  it  is  further  assumed  that  the  failed  unit  is  replaced  by  a  new  unit. 
To  satisfy  the  relibility  requirement,  the  set  of  distribution  functions  of  time- 
to -fir  St -failure  and  time-between-failures  must  be  found.  For  example,  if  an 
80%  confidence  limit  is  specified,  then  the  lower  20%  points  of  all  failure  time 
distributions  must  exceed  the  required  value.  (See  Figure  II-ll,  ) 


Figure  11-10  Fl0«»  >1-11 


To  illustrate  we  will  consider  a  two-unit  system  with  given  failure  time 
distributions  F^(t)  and  F^lt).  See  Figure  11-12. 


time  *1  »2 


n-t 


FIgura  11-12 


To  generate  a  eample  of  time-to-failure,  one  eelects  a  eample  value  (p) 
from  a  tmiform  distribution  betweeon  (0,  1)  and  determines  the  appropriate  failure 
times  tj  andt^s*  shown  above.  System  tlme-to-first>fallure  is  the  smaller  of  the 
two  values  t^,  .  For  this  example  12^^!*  time  the  second  unit  is 

replaced  by  a  new  unit.  The  failure  distribution  curve  for  Unit  #2  is  now  dis¬ 
placed  by  an  amount  t2  •  as  shown  in  Figure  11-13. 


e 
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(»)« 


♦  >»* 


r*(t) 


Fj  (») 


The  failure  time  distribution  for  Unit  #1  is  the  conditional  probability  that 
the  unit  has  not  failed  in  (0,  t2  )  and  fails  at  t  where  t<t2  .  This  conditional 
distribution  F^^ft)  is  given  by; 


Fi(t) 


Fj(t)  -  Fj{t2) 

-T-rrptp- 


t>t 
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To  obtain  the  time -to- second  failure,  we  again  select  a  deviate  from  a 
uniform  distribution  and  proceed  as  above,  thus  obtaining  a  set  of  time-to- second- 
failure.  The  smaller  of  the  two  is  the  time-to-second-failure  of  the  system. 
Proceeding  in  a  similar  nianner  times-to-third,  -fourth,  -fifth  failures,  etc., 
are  obtained,  and,  accordingly,  the  times-between-failures.  This  process  is 
repeated  until  a  sufficient  number  of  times  to  first  failure,  times  to  second,  etc. 
are  obtained.  Thus,  the  empiric  time-between-failure  distributions  can  be  de¬ 
termined  to  any  desired  degree  of  confidexice.  To  meet  the  requirements,  the 
lower  20%  value  ot  all  moan-time-between-failures  must  exceed  that  value 
specified. 
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APPENDIX  III 

PROBABLE  ERROR  IN  STATISTICAL  SOLUTION  OF  AREA  INTEGRAL 

lU.  1  MONTE  CARLO  MODEL  FOR  DETERMINING  AREA 

Thia  appendix  will  determine  the  probable  error  resulting  from  a  given 
nximber  (n)  of  replications  when  Monte  Carlo  is  used  to  ascertain  an  area. 
Figure  111-1  illustrates  the  area  (shaded)  to  be  determined  by  Monte  Carlo 
techniques. 


Figure  IIH 

Let  B  be  the  area  under  f(x)  between  the  limits  a  and  b  ,  thus: 

b 

B  -  y  f(x)dx  (Ul-l) 

a 

The  area  A  under  the  rectangle  is  given  by: 

A  •  (b  .  a)  h  (lU-Z) 


To  determine  the  area  under  f(x)  using  M<mte  Carlo,  mie  proeea^  as 
Icdiosts: 
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From  a  vinlform  distribution  select  a  value  Xj  on  the  interval  from  a  to 
b  .  In  a  similar  manner  pick  from  the  uniform  Retribution  on  the  interval 
frmn  0  to  h  .  These  two  values  define  a  point  Pj  with  coordinates  Pje(X^,  Y|). 
This  process  is  continued  n  times  until  one  obtains  Pj^  •  (X^,  Yj^).  The  limit 
of  the  ratio  of  the  number  of  pointe  falling  in  B  to  the  total  number  of  points 
(n^)  defines  the  percentage  probability  as  n^-*«  or: 

P  ■  llpi~  {UI-3) 

"t 


n-« 

To  obtain  the  area  B»  one  eimply  takes  the  product  PA.  If  the  number  of 
pointe  is  infinite.  PA  ■  B  ,  the  exact  area  f(x).  However,  since  n^  is  finite,  a 
certain  error  is  introduced.  To  find  this  probable  error  as  a  function  of  n*  one 
proceeds  as  follows. 

Define: 

P  a  probability  of  falling  in  B 
1  •  P  ■  probability  of  not  falling  in  B 

The  probability  of  obtaining  r  successes  in  n  tries  can  be  obtained 
from  the  binomial  distribution: 


P(n,r)  ■  cj'p'd  -  P)""' 


(m-4) 


or 


The  average  number  of  successes,  r,  is  given  by: 


P 


P  -  H 


(m-5) 


For  the  binomial  distribution  the  standard  deviation  a  is  defined  by 


a  •  <^nP(l  -  P) 


(m.7y 


$■* 


It  follows  that  one  standard  deviation  in  is; 

n 


and  the  probable  error  (percent)  in  P(n,  r)  is 


(in-8) 


(m-9) 


As  seen  from  equation  (III-9)  the  larger  P  and  n,  the  smaller  the 
probable  error.  Therefore,  the  size  of  the  enclosing  reference  area  should  be 
kept  as  small  as  possible. 


APPENDIX  IV 


DERIVATION  OF  THE  FREQUENCY  FUNCTION  OF  A  FUNCTION 
OF  TWO  INDEPENDENT  RANDOM  VARIABLES 

IV.  1  INTRODUCTION 

This  appendix  considers  the  mathematics  of  a  problem  which  arises  in 
reliability  theory,  viz. ,  the  determination  of  the  distribution  of  the  parameters 
of  a  system  given  the  distributions  of  the  parameters  of  the  components  that 
make  up  the  system.  That  is,  it  is  the  problem  of  determining  the  distribution  of 
a  function  of  several  random  variables  of  known  distributions. 

IV.  2  DERIVATIONS 

Given  the  two  random  variables  x  and  y  ,  their  joint  frequency  function, 

p  (4  fj)  ,  and  some  fimction  relating  a  variable  z  and  the  two  variables  x  and 
X|  y  f 

y,  viz., 

Z  =  f(x,y)  (IV- 1) 

then  z  will  also  be  a  random  variable  whose  distribution  is  determined  by  the 
distr^ution  of  x  and  y.  We  will  develop  methods  for  determining  the  fre¬ 
quency  function  of  the  distribution  of  z. 

The  distribution  function  of  z,  F  (L)  ,  will  simply  be  the  integral  of 

P  (Q  over  that  region  of  the  {  -  ^  plane  for  which 
Xf  y 

1(4,1;) <  ;  (iv-2i 

That  is, 

and  the  frequency  function  of  z,  p,(L),  will  be  the  derivative  of  F  (()  with 
respect  to  ;  i.  e. , 

^Subscripts  are  used  throughout  this  memoramlum  to  identify  fuaeUrnks;  tihe 
subscripts  are  not  to  be  regarded  as  variables. 


lV-1 


p.(t) 


av-4) 


We  will  aaeume  that  equation  (IV- 1)  and  the  Inequality  (lV-2)  may  be 
solved  explicitly  for  y  and  ^  respectively,  and  that  the  solution  of  equation 
(IV- 1)  for  y  is  single-valued  and  differentiable;  thus 


and 


f(x,y)  =  z 

(a) 

r^(x,  z)  ■  y 

(b) 

f(§,9)<  ; 

(a) 

when 

(b) 

} 

} 


(IV-5) 


(IV-6) 


The  inequality  (IV-6b),  corresponding  to  the  inequality  (IV-6a),  may 
have  either  sense,  as  noted,  or  it  may  have  different  senses  in  different  regions 
of  X. 


Then  Equation  (IV- 3)  may  be  written; 

+•0  f’^(|,t) 

-«0 


2-  “  7  >  ,  V  when  f(4,  j^)<i 


(IV-7) 


“  i*  I,  (IV-8) 

The  notatiM  £“  is  the  inversion  of  £(x,  y)  with  respect  to  y;  it  should  he 
confused  with  a  partial  derivative  with  respect  to  y.  • 


JV-2 


If  X  and  y  are  statiatically  independent,  which  will  be  aeaumed  through¬ 
out  the  balance  of  thie  appendix,  the  Joint  frequency  function  of  x  and  y  will 
be  simply  the  product  of  the  frequency  functions  of  x  and  y: 

Px,y(^.7)  Py<^> 

Throughout  the  balance  of  this  memorandum  the  following  assumptions 
(some  of  which  are  implicit  in  the  preceding)  will  be  made: 

1.  PJO  and  p  (m)  are  continuous. 

X  y ' 

2.  The  function,  f(x,  y),  and  its  inversion  with  respect  to  y,  f~  ^(x,  z), 
are  continuous,  single -valued,  and  differentiable. 

3.  X  and  y  are  statistically  independent. 

4.  X,  y,  and  z  are  real. 

When  X  and  y  are  statistically  independent,  equations  (IV-7)  and  (lV-8) 
reduce  to 


1.  When  (ftt)  when  f(^,^)<  t 


-t-oe 


«■  «.t) 


F.d)  =  J  y  Py(y)  dy  d« 


+  W 


j*  p^(i)  Fy[f‘^<,C)ldf 


r-1 


2.  When  ^>£’(1.1)  when 


+  •0  +«0 

*■.<»  ■  j”  >’xt>  J 


(IV- 10) 


(IV- 11 


(IV- 12) 


but  for  genoral  form*  of  the  frequency  functiona,  *ad  Py(|^). 


1.  If 

*  ■  X  +  y 

(IV-20) 

y  «  z  -  X 

(IV-21) 

Then 

when 

z<  t  1 

y  <  c  *  *  J 

(IV-22) 

Therefore 

+•0 

P^U)  ■  j‘p^(«)  Py(;  -  K) 

•  «» 

(IV-23) 

which  is  the  familiar  convolution  integral 

%(t)  •  J  Px(^)  PyU  - 

(IV-24 

2.  If 

z  ■  X  •  y 

(IV-25) 

y  .  £ 

(IV-26) 

-  -  Then . 

. . .  ,  *  <  t 

when 

y<W.  *>0  1 

y  >  1/*  x<  0  J 

**  t 

(IV-27) 

Therefore 

0 

p^d)  ■  -  Jpx<^)Py  k  h  i  ^ 

♦  K‘<'  'r  i  ft  i 

(lV-28) 
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-  y  P',j(f)  {l-Fylf-^(^,t)l}  « 

•  «0 

(IV- 13) 

=  1-3  P3,(4)Fy[r^($.t)l  di5 

(IV- 14)" 

t 

Equation  (IV-4)  for  the  frequency  function  of  becomes,  for  the 

two  cases: 


1.  When  yj  when  f(|,^)<t 

•  «e 

=  3  P^«)  (IV.16) 

(IV- 17) 

(IV- 18) 

We  will  now  consider  a  number  of  specific  cases,  that  is,  a  number  oi 
specific  forms  of  the  function 

*»f(x,y)  (IV- 19) 

IV-5 


2.  When  7>  4)  when 


+  •0  ’  +to 

p.(4l  =  ^  y  P,«)  ^  Py(?) 


+  •0 


r  1  8f*  (£.0 


/  — 


r 


P.U) 


0  -«■« 


•  'y  *  ^x'^'-yUl-r 

.M  0 

<IV-29) 

which  reduce*  to  the  following  form,  which  hae  been  called  the  logarithmic 

convolution: 

m 

“  J  Py  (t)  ^ 

««0 

(IV- 30) 

3.  If 

(IV.31) 

y  ■  zx 

(lV-32) 

Then 

*  <  t 

when 

y  1  ;x  X  >0  1 

y  >  t  X  X  <  0  j 

(IV.33) 

Therefore 

0 

p,.(;)  *  -  J  P,(^)  PyU$)-|^ 

•  «0 

+«» 

(IV-34) 

+  J  P,(^)  Py(U) 

0  « 

A  ^ 

•  jpx<^>  Py(U)^C)‘*^  +  VPx<5>  Py<tl)C^ 

•m  Q 

(IV- 35) 

* 
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m 

(IV.36) 
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IV.6 
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4.  If  z  ■  (IV- 37) 

(Obviously,  z  >  0  if  x  and  y  are  real;  the  latter  assumption  is 
made  throughout  this  memorandum.) 


Then 

when 

Therefore,  when  z  >  0 


.  In  X  *  X  +  y 
y  =  In  z  -  X 


z  <  4 
y  <  ln4  -  X 


(IV- 38) 
(IV- 39) 
(IV-40) 


Pjt)  »  y  )  Py(ln  t  -  4” )  (In  4  -  4 )  d4 


^  j'p^(4)  Py(l«t  - 


•  •  P-(4)  ®  1 

[p^  *  Pyl  (In  4) 


4<  0  1 

t>  0  J 


(IV-41) 

(IV- 42) 

(IV-43) 


where,  in  Equation  (IV-43),  the  usual  notation  for  the  simple  convolution  of  the 
frequency  functions  p^  and  p^  is  used. 

5.  If  z  =  x^  x>  0 

(Obviously,  z  >  0  if  x  and  y  are  real. ) 


InB*  y  In  X 


(IV-45) 


y 


In  z 
In  X 


(IV-46) 


IV-f 


Then 

*  <  t 

when 

X  >  1 

y^^ 

0  <  X  <  1 

(IV-47) 


Therefore,  when  (  0 


'‘i.'i*  *  ■  y  Px'«'  Py  (iSi)  n  (t^) 
0 


+  3  Px<5)Py  (t^)  il  (fe 


(lV-48) 


0 

i  (t^)  T51 


P,(0  = 


t<  0 


I  ^ **y  (^) 


(IV -49) 


(IV- 50) 


We  may  formalize  several  generalities  which  are  suggested  by  the 
preceding: 

1.  Heuristically  we  may  condense  the  two  cases  represented  by  equations 
(IV- 16)  and  (IV- 18)  to  a  single  equation: 


p,(t>  -  j‘p,(T)Py[rV.ol 


»CV.t) 


(IV-51) 


rv-8 


Equation  (lV-51)  may  be  regarded  ae  the  generallaed  convolution  of  p  and  p 

X  y 

with  reipect  to  the  function  f;  this  operation  we  will  denote  by  the  operator 

This  convolution  of  p  and  p  then  becomes  an  operator  on  (  ,  and  equation 

X  y 

(IV-51)  may  be  written,  in  operational  form: 

P^iO  =  IPx*fPyl  (t)  <IV.52) 


2.  Consider  a  functional  relation  between  x,  y,  and  z  of  the  form 

z  =  h[f(x,y)]  (IV-53) 

where  x  and  y  are  random  variables  with  frequency  functions  p  and  p  . 

X  y 

We  may  define  another  random  variable,  w  : 

w  =  f(x,y)  (IV- 54) 

whose  frequency  function,  P^(w  )•  will  be 

P^(»)  =  [Px*fPyl  (IV-55) 

Now 

z  =  h(w)  (lV-56) 

which,  we  will  assume,  may  be  solved  for  w  : 

w  =  h’^z)  (IV-57 

where  h~^(z)  is  a  single-valued  function  of  z  , 


The  frequency  function  of  z  in  terms  of  p^  will  then  be 

P,(t)  »  [h"^(;)] 


(IV-58) 


IV-f 


I  I  I  I  alirti 


Therefore 

Pz't)  *  *f  Pyl 

Example  4  above  is  of  this  type  (see  equation  (lV-48)  1 


V 


{IV-59) 


3.  Consider  a  functional  relation  between  x,  y,  and  z  of  the  form 

*  =  f[g(x).  h(y)]  (IV- 60) 


where,  again,  x  and  y  are  random  variables  with  the  frequency  functions  p 
and  Py..  Define  the  two  new  random  variables; 


u  =  g(x) 

V  =  h(y) 


(IV-61) 


and  ass\une  that  equations  (IV-61)  can  be  inverted  to  give  the  single -valued 
functions  of  u  and  v  ; 


X  =  g"^u) 
y  =  h"^(v) 


(lV-62) 


Then  the  frequency  functions  of  u  and  v  will  be 

=  Pxt8~^^)l‘  ■  1 


(IV-63) 


Substituting  equations  (IV-61)  into  equation  (IV-60)  gives 


s  »  f(UtV) 


(IV-64) 


IV-ifl 


•o  that 


P,(t)  »  [  Pu  *£P^1  (t) 


(IV-72) 


Therefore 


P^U)  =  {{Pjt*gPyl  *fP^)(t) 


(IV-73) 


5.  Again,  let  x,  y,  u,  and  v  be  random  variables  with  frequency  func- 

t 

a  =  f  I  g  (x.  y).  Mu,  v)  ] 


tions  ,  Py  ,  p^ ,  p^  ,  and  let 


If  we  let 


•  =  g(x.y) 

t  =  h(u,v) 


} 


then  s  and  t  are  random  variables  with  frequency  functions 


)  »  [  P^  ♦g  Py  1  (d )  1 

Pt(T)  =  [p^\Pyl  (T)  J 


X  is  a  function  of  s  and  t  : 

8  s  f  (s.t) 

and  is  a  random  variable  with  the  frequency  function 

VgiO  ■  fP*  *|PtJ 

Therefore 

p,.(t)  ■  {  IPx  *g  Pyl  *£  IPy  ♦h  PyJ  ^ 


{IV-74) 


(IV-75) 


{IV-76) 


(IV-77) 


(IV-7i) 


(IV-791 


IV-12 


IV.  3  SUMMARY 


If  we  regard  p-^  a«  an  operator,  i.  e. ,  if  we  regard  result 

of  performing  the  operation  p  on  ,  we  may  stunmarize  the  results  of  the 

z 

preceding  in  operational  notation. 

Let  X  ,  y  ,  u  ,  V  ,  be  random  variables  with  frequency  fxmctions  defined  by 
the  operators  p^  ,  p^,  p^  ,  p^  ,  and  let  z  be  a  function  of  the  variables  x  ,  y  , 
u  ,  V  ,  viz. , 


z  =  f(x,  y,  u,  v) 


(IV-80) 


Then  z  will  be  a  random  variable  whose  frequency  fiinction  will  be  defined  by  the 

operator  p  .  We  will  tabulate  the  form  of  the  operator  p  for  various  forms 
z  z 

of  the  function,  f(x,  y,  u,  v)  (i.  e. ,  the  operator  f ). 


First,  we  will  define  the  generalized  convolution  operator  with  respect  to 
the  function  f  of  two  variables,  viz. ,  ,  as  : 

«o 


aCV,;) 


lPx*£PyH4)=  r  ^ -  dT 

•  «d 


,-l, 


where  f^  (x,  z)  is  the  inversion  of  f(x,  y)  with  respect  to  y;  i.  e.  , 


given: 

then: 

The  tabulation  is  as  follows 


2  =  f(x,y) 


y  =  fy^x.z) 


(IV-81) 


{IV-82) 


1.  If 

z  =  f(x,y) 

{IV-83) 

Pz  =tPx*fPyl 

(IV-84) 

2.  If 

N 

II 

& 

Ms 

X 

(IV-85) 

IV-13 


4 

=  [DX[p^*^Pyl  (h"^)]  (IV.86) 

where  D  la  the  differential  operator:  i.  e. ,  ' 

D  [4>  (x)]  3  4x)  a  ^  4>(x)  (IV- 87) 

The  cross  in  equations  (lV-86)  indicates  multiplication  in  the  following  sense; 

Given  two  operators,  and  L2  ,  and  the  variable  (or  function,  or  operator), 

z  ,  then 


Lj  L2{z)s  Lj[L2(z)] 

(IV.88) 

Lj  X  L2(z)4Li{*)1  X  [L2(z)] 

(IV-89) 

That  is,  the  notation  of  equation  (IV-89)  indicates  the  successive  application  of 
the  operators  and  :  first  the  operation  is  performed  on  z  ,  and 

then  the  operation  Lj  is  performed  on  the  result.  The  notation  of  equation 
(lV-89)  indicates  the  product  of  the  results  of  performing  each  of  the  operations 
and  L2  on  z  . 

3.  If 

z  =  f  [  g(x),  h(y)] 

(IV-90) 

=  {[  (g'^)'Pj£(g'Sl*£  [  (h'^)'py(h‘^  ]•} 

=  [{[DXp^l  (g-^}*j  {[DXpy]  (h-^)}] 

►  (IV.91) 

4.  If 

z  =  f[g(x,  y),  u] 

(IV-92) 

Pz  =  ttPx*gPyl  *fPu> 

(IV.93) 

5.  If 

z  =  f[g(x,  y),  h(u,  v)l 

(IV.94) 

P*  =  {tPx  *gPyl  *£  tPu  \Pj) 

(IV.95) 

IV-14 


In  the  following  particular  cates,  the  convolution  takes  the  forms  shown: 


6. 

If 

z  a  f(x,y) 

a  X  +  y 

go 

(IV-96) 

Pj(C)  =  [Pjj^fPyl  (t)  = 

j’px<^)PyU  -  ■") 

(IV-97) 

7. 

If 

*  =  £(x.y) 

=  X  •  y 

(IV-98) 

PgU)  =  [P3j*fPyl  (0  = 

fpx<">  Py  (4) 

(IV-99) 

8. 

If 

Z  a  f(x,  y) 

s  1 

X 

go 

(IV- 100) 

Pa(4)  =  tP^  *f  Py]  (t)  = 

^PX^^^  Py(l  T)|t|  dT 

-M 

(IV-101) 

9. 

If 

*  =  f(x.y) 

(IV- 102) 

10. 


PjtO  =  tPx  *f 


(IV- 103) 


Z  a  f(x,  y)  a 

P,(t)  «  lPx*fPyl  U) 
r  0 


*>0  (IV- 104) 


-fe- 

^  A 


i<  0^ 

no 


(IV- 105) 


IV-15 


Th«  preceding  it  predicated  on  the  following  attumptipne: 

1*  P  •  P  •  P  •  P  &re  continuous. 

•"x  •  r-y  •  ^u  •  ry 

2.  f.  g,  h,  and  their  Invertifins,  f'^,  g"^,  h’^,  are  continu 
tingle -valued,  and  differentiable. 

3.  X  ,  y  ,  u  ^  and  v  are  ttatittically  independent. 

4.  X  ,  y  ,  u  ,  V  ,  and  z  are  real. 


